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Abstract 



The vacuum of quantum electrodynamics is unstable against the formation of many- 
body states in the presence of an external electric field, manifesting itself as the 
creation of electron-positron pairs (Schwinger effect). This effect has been a long- 
standing but still unobserved prediction as the generation of the required field 
strengths has not been feasible so far. However, due to the advent of a new gener- 
ation of high-intensity laser systems such as the European XFEL or the Extreme 
Light Infrastructure (ELI), this effect might eventually become observable within 
the next decades. 

Previous investigations of the Schwinger effect led to a good understanding of the 
general mechanisms behind the pair creation process, however, realistic electric fields 
as they might be present in upcoming high-intensity laser experiments have not been 
fully considered yet. Actually, it was only recently that it became possible to study 
the Schwinger effect in realistic electric fields showing both temporal and spatial 
variations owing to the theoretical progress as well as the rapid development of 
computer technology. 

Based on the equal-time Wigner formalism, various aspects of the Schwinger ef- 
fect in such inhomogeneous electric fields are investigated in this thesis. Regarding 
the Schwinger effect in time-dependent electric fields, analytic expressions for the 
equal-time Wigner function in the presence of a static as well as a pulsed electric field 
are derived. Moreover, the pair creation process in the presence of a pulsed electric 
field with sub-cycle structure, which acts as a model for a realistic laser pulse, is 
examined. Finally, an ab initio simulation of the Schwinger effect in a simple space- 
and time-dependent electric field is performed for the first time, allowing for the 
calculation of the time evolution of various observables like the charge density, the 
particle number density or the number of created particles. 
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Chapter 1 

Introduction 



Ever since the publication of 'Philosophiae Naturalis Principia Mathematica' by Sir 
Isaac Newton in 1687, people have started to believe that the world around us is 
governed by the action of forces. Nowadays, more than 300 years later, we are 
rather convinced that there are four fundamental forces in nature: The electromag- 
netic force, the weak force and the strong force are properly formulated in terms of 
quantum field theories whereas the gravitational force is still best described by the 
classical theory of general relativity. 

The force of interest in this thesis is the electromagnetic force, which is well 
known from everyday life: Phenomena ranging from the radiation of light over 
all home and personal electronic equipment to high-tech devices in industry and 
research are based upon electromagnetic radiation and electric currents, most no- 
tably composed of electrons. It is sometimes appropriate to describe these phenom- 
ena by the classical theory of electromagnetism [1], however, a quantum description 
becomes essential as soon as phenomena such as laser emission are considered. 

The quantum theory of electromagnetism is called quantum electrodynamics 
(QED) and has been developed during the first half of the 20tli century [2, 3, 4], 
which finally culminated in the Physics Nobel Price-'^ awarded to J. Schwinger, 
R. Feynman and S. Tomonaga in 1965. Since then, this theory has been tested 
experimentally to a very high precision in a regime which is characterized by high 
energies and low intensities. In this high-energy regime, which is realized in particle 
accelerator experiments, perturbation theory is a well-suited tool in order to cal- 
culate observables such as the anomalous magnetic moment of the electron [5, 6], 
which is one of the most precisely known physical constants nowadays. 

There is yet another regime which is characterized by low energies and high 
intensities and, therefore, experimentally not accessible in accelerator experiments. 
This strong-field regime, however, might become accessible in the near future by 
the advent of a new generation of high-intensity laser systems. Physicists pin their 
hopes on two different technologies: X-ray free electron laser (XFEL) systems such 

^cf. http: //nobelprize . org/nobel_prizes/physics/laureates/1965/ 
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Figure 1.1: Left: Accelerator experiments are characterized by high energies but low inten- 
sities whereas laser-based experiments operate at low energies but high intensities. Right: 
The attainable intensity in optical high-intensity laser systems has continuously increased 
since the advent of the chirped pulse amplification (CPA) technology in 1985 [7]. The 
figures are taken from [8, 9]. 

as the European XFEL^ or the Linac Coherent Light Source LCLS^ on the one hand, 
and optical high-intensity laser systems such as the Extreme Light Infrastructure 
ELI^ or the High Power laser Energy Research facility HiPER^, which are based on 
chirped pulse amplification (CPA) technology [7], on the other hand. 

The non-trivial and non-perturbative structure of the QED vacuum plays a cru- 
cial role in the strong-field regime: Due to the fact that quantum field theories allow 
for virtual particles being produced according to the energy-time uncertainty prin- 
ciple, the vacuum cannot be regarded as totally empty. In fact, the vacuum has to 
be considered as a polarizable medium with refractive and absorptive index. Con- 
sequences are rather minor in the high-energy regime: The vacuum is considered as 
the ground state of the system and observables are measured with reference to this 
state. In the strong-field regime, however, the non-trivial vacuum structure has a 
significant effect: Whilst the direct interaction between photons is forbidden by C/(l) 
gauge symmetry, photons might interact indirectly via vacuum fluctuations. Con- 
sequently, depending on the laser's frequency and intensity, various effects ranging 
from photon-photon scattering over vacuum birefringence to vacuum pair creation 
might arise. For recent reviews on various quantum vacuum experiments see [9, 10]. 

^cf. http://www.xfel.eu/ 

^cf. http: //Icls . slac . stanf ord. edu/ 

^cf. http : //www . extreme- light - infrastructure . eu/ 

'"'cf. http : //www . hiper- laser . org/ 
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In this thesis I focus on vacuum pair creation (Schwinger effect), which has 
been a long-standing but stiU unobserved prediction of QED [11, 12, 13]: The QED 
vacuum becomes unstable against the formation of many-body states in the presence 
of an external electric field, manifesting itself as the creation of electron-positron 
pairs. Subsequently, it turned out that the Schwinger effect shows two characteristic 
features which make its theoretical treatment both challenging and valuable in view 
of current problems in modern physics: On the one hand, the pair creation process is 
non-perturbative in the electromagnetic coupling so that QED perturbation theory 
fails. On the other hand, the system's time evolution is a non- equilibrium process 
and belongs as such to the lest-well understood branches of quantum field theory. 
Accordingly, a careful investigation of the Schwinger effect might entail important 
lessons for related phenomena such as the formation and thermalization of the so- 
called quark gluon plasma or facilitate the search for hidden sector particles [14. 15, 
16]. 

In order to cope with this challenging non-equilibrium and non-perturbative 
problem, various methods have been invented, which might be roughly divided into 
two groups: On the one hand, there are various semi- classical methods such as 
WKB approximation [17, 18, 19, 20] or instanton techniques [21, 22, 23, 24]. These 
methods relate the imaginary part of the effective action for a given background 
field to the vacuum persistence amplitude, which in turn can be related to the 
pair creation rate [25, 26]. On the other hand, there are various quantum kinetic 
m,ethods such as the quantum Vlasov equation [27, 28, 29, 30] or the Wigner function 
formalism ]31, 32, 33], which account for both the pair creation and the subsequent 
transport process. In addition to the pair creation rate, quantum kinetic methods 
yield valuable phase-space information such as the momentum distribution. 

Initially, the Schwinger effect has been considered for static electric fields [13] 
as well as simple time-dependent electric fields [17, 18]. Since then, various inves- 
tigations have been conducted for more complicated time-dependent electric fields 
[34, 35, 36, 37, 38, 39, 40], space-dependent electric fields [22, 24, 25, 41, 42] as well 
as collinear electric and magnetic fields [43, 44, 45], resulting in a good understand- 
ing of the general mechanisms behind the pair creation process by now. However, 
realistic field configurations as they might be present in upcoming high-intensity 
laser experiments have not been fully considered yet. Most notably, there are two 
aspects which have not been properly taken into account so far: First, a realistic 
electric field shows both spatial and temporal dependence. Accordingly, a subtle 
interplay between the different scales might occur, resulting in a modification of 
both the pair creation and/or the transport behavior. Secondly, the magnetic field 
should be taken into consideration in a realistic manner as well. 
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This thesis is devoted to take a further step towards an accurate description of 
the Schwinger effect in upcoming high-intensity laser experiments. To this end. I fo- 
cus on several aspects of the Schwinger effect in space- and time-dependent electric 
fields, most notably the interplay between spatial and temporal inhomogeneities. 
The Wigner function formalism and the closely related quantum Vlasov equation 
are adopted throughout as they provide adequate methods to study the influence of 
space- and time-dependent electric fields on the Schwinger effect. 

The thesis is organized in the following way: In Chapter 2, I give a short overview 
of the Schwinger effect including the current state of research and the objectives for 
this thesis. The theoretical basis for the subsequent investigations is provided in 
Chapter 3 by reviewing both the covariant and the equal-time Wigner formalism. 
In Chapter 4, the Schwinger effect in time-dependent electric fields is discussed in 
some detail. First, I give a short review on the derivation of the quantum Vlasov 
equation and prove its equivalence to the equal-time Wigner formalism in the limit 
of a spatially homogeneous, time-dependent electric field. Afterwards, the analytic 
solution of the quantum Vlasov equation for the constant electric field as well as 
for the pulsed electric field is given. Moreover, the Schwinger effect in a short pulse 
with sub-cycle structure is studied. In Chapter 5, the Schwinger effect in space- and 
time-dependent electric fields is investigated. After introducing a 1 + 1 dimensional 
model of QED, I present some solution strategies for the corresponding equation 
of motion. The remainder of this chapter is then devoted to the discussion of the 
results of an ab initio simulation of the Schwinger effect in a simple space- and 
time-dependent electric field. In the final Chapter 6, some concluding remarks as 
well as an outlook on potential future investigations are given. 

The thesis concludes with two appendices: In Appendix A, the details of diverse 
derivations and calculations are given which have been deferred from the main text 
in order to improve readability. In Appendix B, the finite difference scheme for 
solving hyperbolic PDE systems is presented upon which the numerical simulations 
are based. 



Chapter 2 

Overview and objectives 



This chapter is dedicated to giving an overview of the major developments in the 
field of strong- field QED with a special focus on the Schwinger effect. In Section 2.1. 
I concentrate on the early milestones which have been the basis for many ensuing 
investigations. In the subsequent Section 2.2, I focus on quantum kinetic methods 
and briefly review some main results. In Section 2.3, I draw my attention to some 
current research topics in the field of vacuum pair creation. Finally, I formulate the 
objectives of this thesis in Section 2.4. 

2.1 Schwinger effect 

The history of the Schwinger effect started shortly after the invention of quantum 
mechanics in the mid 1920s. It was P. Dirac who first set up the relativistic wave 
equation for the electron in 1928 [46]. In the first instance, the negative energy 
solutions of the Dirac equation exhibited an interpretational difficulty: Why does 
a positive energy solution not decay into a negative energy solution by continuous 
emission of photons? This puzzle was first solved by the Dirac sea picture, according 
to which the vacuum of the theory consists of a completely filled, negative energy 
band as well as an empty, positive energy band, which are separated by an energy 
gap of 2m. ^ In this picture, a hole in the negative energy band is interpreted as the 
electron's antiparticle, the positron.'^ 

Subsequently, the Dirac equation was solved in the presence of a static electric 
field Eq by F. Sauter in 1931 [11]. He found that the positive and negative energy 
bands are bend by the electric field so that a level crossing occurs. In the Dirac sea 
picture, this means that a negative energy electron is allowed to tunnel through the 
energy gap to the positive energy band, leaving a hole behind. In modern language, 

^Note, that the Dirac sea picture is not apphcable to relativistic theories of scalar particles, even 
though the negative energy problem exists there as well. Consequently, the notion of the Dirac 
sea should be dropped in favor of a quantum field theoretical treatment of the negative energy 
solutions. 

^Natural units ft = c = 1 are used throughout the thesis. 

■'Since the positron has not been observed until 1932, it was the proton which was first inter- 
preted as the electron's antiparticle. 
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this result states that there is a non-vanishing probabihty for electron-positron pairs 
to be created in the presence of an electric field. As the pair creation process might 
be considered as tunneling process, it is exponentially suppressed though; 

P[e+e-]^exp(-^j ' ^^.l) 

with the scale set by the electron mass. Accordingly, electric fields of the order of 
the critical field strength E„ are needed to overcome the exponential supression: 

E„ = — ~ 10^*^ V/m . (2.2) 

e 

These early investigations indicated that the Dirac equation is by no means a 
single-particle equation but should rather be treated in a quantum field theoretical 
framework in order to account for the possibility of particle creation. Accordingly, 
the vacuum should be considered as polarizable medium due to permanent vacuum 
fluctuations. As a consequence, the linear Maxwell's equations have to be replaced 
by more complicated non-linear equations. 

This peculiarity has first been pointed out by W. Heisenberg and H. Euler in 1935 
[12]. They derived an integral representation of the one- loop effective Lagrangian 
'^EH which accounts for the coupling of a static electromagnetic field F^^'^ to the 
electron vacuum loop. The leading order correction to Maxwell's Lagrangian is 
then given by quartic terms in the field strength tensor, resulting in non-linear 
corrections to Maxwell's equations: 

C% = -:F+ ^^^[4^^ + 7g'] + 0{F^) . (2.3) 

Here, J" = iFf^'F,,^ = ^(Bo^ - Eg) and Q = IF^'" F,,^ = -Eq ■ Bq denote the in- 
variants of the electromagnetic field. Since then, the effective Lagrangian became a 
powerful tool in describing non-linear effects such as photon-photon scattering, vac- 
uum birefringence or vacuum pair creation. Note for completeness that an integral 

(2) 

representation of the two-loop effective Lagrangian accounting for a photon 

correction to the electron vacuum loop, has been derived as well [47]. For a recent 
review on the topic of Euler-Heisenberg effective Lagrangians see [48]. 

In his seminal work, J. Schwinger was the firsf* to calculate the imaginary part 
of the one-loop effective Lagrangian in the presence of a static electric field Eq 

''Actually, it had been pointed out by W. Heisenberg and H. Euler that the imaginary part of 
the effective action could be calculated analytically by taking the pole sum in their expression, 
however, they did not give the final result [12]. 
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exactly [13]. Since then, vacuum pair creation by electric fields has been widely 
referred to as the Schwinger effect: 



electron-positron pairs are created in the presence of a static electric field Eq.^ 
This statement, however, is not totally correct as this expression should rather be 
identified with P[vac.], that is the rate at which the vacuum decays. In fact, the 
vacuum decay rate is always larger than the pair creation rate, which is given by 
the first term in this series only [25, 26]: 



As a matter of fact, the difference between Eq. (2.4) and Eq. (2.5) is rather minor in 
the sub-critical field strength regime whereas it becomes important for super-critical 
fields. Note for completeness that Eq. (2.5) can be generalized for arbitrary spins 
and space dimensions [49]. 

Subsequently, there has been further research on the Schwinger effect in time- 
dependent electric fields E(f) as well. These investigations have been largely moti- 
vated by studies on atomic ionization in either static [50] or sinusoidal [51] electric 
fields. There is indeed a close analogy in spirit between vacuum pair creation and 
atomic ionization: Most notably, both of them can be described as tunneling process 
which is facilitated by an applied electric field. ^ Accordingly, methods and notions 
which were used in the theoretical treatment of atomic ionization have also been 
applied in the investigation of the Schwinger effect. 

Atomic ionization in a sinusoidal electric field E(t) = Eq sin(aji) is one prominent 
example: This is, in fact, a multiple time scale problem, with the relevant scales 
given by the frequency uj as well as the tunneling frequency ujt, which is related to 
the characteristic time of a tunneling event. The Keldysh adiabaticity parameter 7 

^Eveii J. Schwinger referred to this expression as the pair creation rate [13]. 
Amongst others, this can be seen from the expressions for the pair creation and ionization rate, 
respectively. The main difference is the scale, which is orders of magnitudes smaller for atomic 
ionization, whereas the general structure is the same {£b denotes the binding energy): 





This expression has often been identified with J\f[e'^e ], that is the rate at which 






and 
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then determines the dominant scale of the problem: 



leEol ^ ' 

with Eb denoting the binding energy. For 7^1, the problem is dominated by ojt so 
that the electric field might be well considered in an instantaneous approximation. 
As the ionization is due to a tunneling event, this type has been termed instanta- 
neous or non-perturbative ionization. For 7^1. however, the temporal variation 
of the electric field is too quick so that no tunneling event is possible. In this case, 
the ionization is due to the absorption of a certain number of photons so that this 
effect has been named multi-photon ionization. 

The notion of the Keldysh parameter has been adopted in strong-field QED as 
well. Again, 7 measures the ratio between the characteristic frequency of the electric 
field uj and the characteristic tunneling frequency ojt'- 

7=- = ^- (2.7) 
OJT leEol ^ ' 

As for atomic ionization, the non-perturbative regime is characteristed by 7 <C 1 
whereas 7^1 denotes the multi-photon regime.^ As a matter of fact, an analytic 
expression for the vacuum decay rate has been derived which interpolates between 
those two regimes [17, 19]. Considering the pair creation rate which is in the sub- 
critical field strength regime again very close to the vacuum decay rate, one obtains: 



exp — 7 4C 1 (instantaneous) , 

^[^^^~]^{ n^.X/J (2-8) 
2^) 7 » 1 (dynamical) . 

This result nicely illustrates the meaning of the Keldysh parameter: On the one 
hand, the characteristic exponential supression of a tunneling event is found for 7 <^ 
1. On the other hand, the result for 7^1 corresponds to n-th order perturbation 
theory, with n being the minmum number of photons to be absorbed in order to 
overcome the threshold energy for pair creation nuj > 2m. 

Moreover, the pulsed electric field E(i) = Eosech^(|:), with t being a measure 
for the pulse length, has been of special interest as it describes an electric field which 
is only switched on for a finite time period [18]. Again, an analytic expression for 
the vacuum decay rate in terms of an integral representation has been derived for 
this type of electric field. I will return to this issue in Section 4.3. 



''in the context of pair creation, this effect is also termed dynamical pair creation. 
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Figure 2.1: Schematic representation of a heavy ion coUision. From left to right: Ap- 
proach of two nuclei; heavy ion collision; formation of the QGP; expansion of the QGP and 
hadronization; hadronic freeze out. The figure is taken from http://www.phy.duke.edu/ 
research/NPTheory/QGP/transport/ 

2.2 Quantum kinetic methods 

The Schwinger effect attracted much attention not only in the field of strong-field 
QED but also in other research areas. One important example is the theoretical 
description of the pre-equilibrium evolution of the so-called quark giuon plasma 
(QGP) in ultrarelativistic heavy ion collision at the Relativistic Heavy Ion Collider 
RHIC^ or at the Large Hadron Collider LHC^. In order to describe the formation of 
the QGP, one prominent model considers the colliding nuclei as passing through each 
other, creating a chromoelectric flux tube which breaks via the production of quark- 
antiquark pairs [52]. Even though this phenomenologically motivated description 
bears a number of shortcomings, most notably it is by no means clear that a non- 
Abelian chromoelectric field can be treated on the same ground as an Abelian electric 
field [53], much effort has been made to find a quantum kinetic description of non- 
equilibrium systems showing the phenomenon of particle creation. 

2.2.1 Quantum Vlasov equation 

Kinetic theory has been very successful in describing transport phenomena in non- 
relativistic theories such as fluid dynamics as well as in relativistic applications 
such as plasma physics. Even the generalization to quantum systems can be done 
consistently. For a review on kinetic theory see [54]. 

The kinetic description of non-equilibrium systems is based on the Boltzmann 
equation^ which describes the time evolution of the one-particle distribution function 
.7-'(x, p, f). In general, this equation reads: 



J^(x,p,0 = [|+x- Vx + p- Vp] ^(x,p,f) =C[^]. 




*cf. http://www.bnl.gov/rhic/ 
cf. http: //public .web . cern. ch/public/ 
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Here, x denotes the spatial variable and p is the momentum variable. Moreover, the 
dot represents a total time derivative and C[J-] denotes the collision term, which is 
unspecified so far. This equation states that the temporal change of the one-particle 
distribution function has three reasons: The particle's motion, the action of forces 
as well as collisions between particles, hi order to include pair creation into the 
kinetic description, one has to account for this by an additional source term <S[jF]. 

The importance of collisions strongly depends on the problem: On the one hand, 
the density of particles is expected to be quite high in heavy ion collision so that col- 
lisions surely have to be taken into account. There have been various investigations 
of this effect by means of relaxation time approximations [55, 56, 57] as well as more 
sophisticated models [58, 59, 60]. On the other hand, the density of particles is ex- 
pected to be rather low when considering the Schwinger effect in the sub-critical field 
strength regime so that collisions might be neglected to a good approximation [61]. 
Taking this into account, the collisionless approximation of the Boltzmann equation 
(quantum Vlasov equation) for the case of a spatially homogeneous, time-dependent 
electric field E(i) reads: 



A first ansatz for S[J-] was based on the Schwinger formula Eq. (2.5) and assumed 
the electron-positron pairs to be created at rest ]28, 52, 62, 63]: 



tum in the orthogonal direction. Note that S[J-] does not depend on J-{p,t) in 
this ansatz. It is clear that the validity of this approximation is restricted: The 
Schwinger term had been derived for a static electric field Eq whereas this ansatz 
assumes a time-dependent electric field E(t). Consequently, the applicability in the 
case of a rapidly varying electric field is not well justified. 

It was only in the mid 1990s that the nowadays widely used source term for the 
Schwinger effect in time-dependent electric fields E(f) was derived. It was first 
shown for scalar QED (sQED) [29] and later also for QED [30] that S[J^] can be 
found from first principles in mean field approxim,ation. In QED, this source term 
reads: 




(2.10) 




(2.11) 



Here, py denotes the momentum in the field direction whereas p_L is the momen- 



t 




(2.12) 



— oo 
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The detailed derivation of this source term with all quantities properly defined is 
given in Section 4.1. Note that this source term indeed depends on J-{p, t) as distin- 
guished from Eq. (2.11). Accordingly, it was pointed out that the Schwinger eff'ect 
possesses various characteristic features: Most notably, the pair creation process 
shows non-Markovian behavior because of memory effects and non-locality in time 
[61, 64]. Additionally, it was shown that the quantum Vlasov equation exhibits the 
features of quantum statistics, i. e. Bose enhancement for scalar particles as well as 
Pauli-blocking for spinor particles [30]. For computational reasons, it turned out to 
be of advantage to rewrite this integro-differential equation as a first order, ordinary 
differential equation (ODE) system [65]: 



Again, the detailed derivation of this ODE system is given in Section 4.1. 

The rather clear and simple structure of the ODE system made detailed analysis 
of the Schwinger effect in any kind of time-dependent electric field E(i) feasible. Ini- 
tially, these investigations dealt with rather simple field configurations such as static 
[64, 65], sinusoidal [66, 67, 68] or pulsed [57, 69, 70] electric fields. Nonetheless, these 
studies revealed remarkable features such as a non-trivial momentum dependence of 
the pair creation process, the phenomenon of pair creation and annihilation in tune 
with the frequency of the sinusoidal field or the strong enhancement of the Schwinger 
effect in electric fields with a temporal variation of the order of the Compton time 
Tc- Additionally, it was shown that the density of created particles is rather low 
in the sub-critical field strength regime so that the omission of the collision term 
C[J-'] seems well justified. Very recently, further investigations of more realistic field 
configurations such as pulsed electric fields with sub-cycle structure have been car- 
ried out [36, 39]. It was shown that the momentum distribution of created particles 
strongly depends on the detailed form of the electric field. 

A further problem which has been studied is the influence of the self-induced 
electric field due to the pair creation process: The applied electric field Eext(i) 
creates charged particles which in turn produce an opposing internal electric field 
Eint(i)- Accordingly, one has to consider the Schwinger effect in the resulting electric 
field E(t) = Eext(i) + Eint(i) in a self-consistent way. Again, this phenomenon has 
been investigated in mean field approximation, where an analytic expression for 
the internal electric field can be derived [27, 65]. Loosely speaking, it consists of 
three contributions: One contribution because of the motion of existing particles 



^(p,i) 



Qip,t)g{p,t) , 

Q{p,t)[l - J'ip^t)] - 2uj{p,t)n{p,t) 
2u;{p,t)g{p,t) . 



(2.13) 
(2.14) 
(2.15) 
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(conduction current), another one due to the continued pair creation (polarization 
current) and a third one because of the necessity of charge renormalization. As for 
the influence of collisions, it has been shown that this backreaction mechanism can 
be well ignored in the sub-critical field strength regime [65, 71]. 

2.2.2 Wigner function formalism 

The phase space formulation of quantum mechanics has first been introduced in 
order to account for quantum corrections to the Boltzmann equation. It turns 
out that this approach is by no means unique but allows for a large number of 
different realizations in terms of, most prominently, the Wigner function [72], the 
Husimi function [73] or the Glauber-Sudarshan function [74, 75]. For a thorough 
introduction see [76]. 

I will focus on the Wigner function W(x, p, t) here, which has been widely 
adopted in the treatment of non-relativistic quantum physics. The phase space 
formalism has the advantage of offering an approach which uses classical language 
in the description of quantum systems as well as allowing direct access to real-time 
dynamics. Even though the Wigner function cannot be considered as a classical 
distribution function, it is in fact an observable which has been measured in various 
optical and atomic setups [77, 78, 79]. 

In a non-relativistic context, the Wigner function is defined as the Fourier trans- 
form of the density matrix. "'^'^ Accordingly, the relative coordinate y is traded for 
the momentum variable p: 



with dF = d'^x d^p/(27r)^ being the phase space volume element. As already men- 
tioned, the Wigner function cannot be considered as a distribution function but 
should rather be termed a quasi- distribution function since it can take negative val- 
ues W(x, p,i) < 0. It is especially this peculiarity which makes its interpretation 
quite challenging. Nevertheless, it still possesses various desired properties, most 

^"For a pure state l^"), the density matrix reduces to p = and the Wigner transform reads 




(2.16) 



The expectation value of any quantum mechanical observable O is then calculated 
according to; 
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notably, it yields the correct quantum-mechanical marginal distributions: 




and 



/ (|^mx,P,i) = (p|p|P>- (2.18) 



The need for an appropriate description of transport processes in astrophysics as 
well as plasma physics applications led to the generalization of the non-relativistic 
Wigner function to a relativistic context [80]. The description of both Abelian gauge 
theories like QED and non-Abelian gauge theories like quantum cliromodynamics 
(QCD), however, entails still another complication; Most quantities such as corre- 
lation functions depend on the chosen gauge whereas all observable quantities must 
not. Due to the fact that the Wigner function should at least in the classical limit 
coincide with the corresponding classical distribution function, which in fact is an 
observable quantity, it is necessary to define the covariant Wigner function in a 
gauge invariant way [81, 82, 83]. 

As my focus lies on the Scliwinger effect, I restrict myself to QED in the following. 
The covariant Wigner operator VV('^)(x, p) is then defined in a similar manner as it's 
non-relativistic counterpart^^ : 



with ^'(x), ^'(x) and A^^{x) being considered as Heisenberg operators. It has to 
be emphasized that the path-ordered Wilson line integral U{A;x,y), which ensures 
gauge invariance. is by no means unique. A physical sensible interpretation of p 
as kinetic momentum, however, forces the integration path to be chosen along the 
straight line [84]: 



the corresponding Wigner operator. Due to the fact that I focus on vacuum pair 
creation, the ensemble average can be replaced by the vacuum expectation value: 



where |Q) denotes the Heisenberg vacuum. Note that the resulting equation of mo- 
tion for W'-'*-' (x, p) forms a tower of coupled equations due to the operator nature 

^""^Note that depending on the situation four-vectors are denoted either by sans serif variables or 
by explicitly writing the Lorentz indices, for example p = p** = (p", p) 




(2.19) 




W(^)(x,p)^(0|w(4)(x,p)|0), 



(2.21) 
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of Af^{x). This infinite BBGKY hierarchy, named after N. Bogoliubov, M. Born, 
H. Green, G. Kirkwood and J. Yvon, has to be truncated for practical reasons. 
Treating the vector potential as classical background, this truncation of the BBGKY 
hierarchy happens automatically at the one-body level. Note that this formalism 
is put on the same footing as the quantum Vlasov equation once all photon con- 
tributions beyond mean field approximation are neglected. Further details on the 
covariant Wigner formalism can be found in Section 3.1. 

The covariant Wigner formalism provides a clear theoretical description, however, 
it is rather inconvenient for various reasons. Consider the following situation in an 
upcoming high-intensity laser experiment: Electron-positron pairs are created in 
the focus of colliding laser pulses after the system had been prepared in the vacuum 
state Subsequently, the created particles move apart from each other with the 
positrons being detected by surrounding detectors. 

This idealized experimental situation has implications on the theoretical de- 
scription: First, in order to trace the system's time evolution, the problem should 
be posed as an initial value problem. Consequently, it is necessary to choose a dis- 
tinguished time t = x*^ and give up explicit Lorentz covariance. Secondly, it is not 
required to stick to a gauge invariant approach once the laser pulse is described as a 
classical electromagnetic field. Accordingly, it is convenient to choose a gauge fixed 
expression for the vector potential A'^(x). 

There have been two approaches to derive such an equal-time Wigner formalism: 
On the one hand, there have been attempts to define the equal-time Wigner function 
W*'^^(x, p,t) as an energy average of the covariant Wigner function [85, 86, 87]: 



W(3)(x,p,i)= / ^W(^)(x,p). (2.22) 



On the other hand, the equation of motion for the equal-time Wigner function has 
also been derived in a non-covariant and gauge-fixed manner from the very beginning 
[31]. In fact, both approaches yield the same equation of motion: 

DtW^^^ = -iD' [tV,^^^^] -im [7°,W(3)] - jiV, W^^^ } , (2.23) 

where Dt, D and 11 denote pseudo-differential operators to be defined later. Further 
details on the equal-time Wigner formalism can be found in Section 3.2. To be 
complete note that this formalism has been applied to sQED as well [88, 89]. 



^^Note that there are even other proposals to detect the Schwinger effect, for instance in photon- 
nucleus colhsions. 
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There have been various investigations on the equal-time Wigner formahsm since 
then: First, it has been shown that this formahsm gives a transport equation of the 
form of Eq. (2.9) in its classical limit [90, 91, 92, 93, 94]. Further investigations 
included loop corrections to the classical limit and showed the necessity of renor- 
malization at higher orders [31, 95]. Additionally, there have been first attempts 
to tackle the problem of pair creation within this formalism [31, 61, 89, 96]. These 
investigations were closely related to studies of the quantum Vlasov equation with a 
source term Eq. (2.11). The close connection of the equal-time Wigner formalism to 
the non-Markovian source term Eq. (2.12) has been pointed out only recently [33]. 

In addition, there have been more general investigations of strong-field effects 
by means of the Wigner formalism in the light-front approach lately. These studies 
showed that the Wigner function incorporates all characteristic features of strong- 
field QED, most notably the mass shift of a particle in a plane wave background as 
well as effects due to multi-photon absorption and emission [97]. 

2.3 Current research topics 

After giving this brief overview of the theoretical description of the Schwinger effect 
by means of quantum kinetic methods, I now turn to some topics which are of special 
interest in the current research. 

2.3.1 Momentum distribution in realistic laser pulses 

The upcoming high-intensity laser facilities will be capable of doing fundamental 
physics in the strong-field regime [9, 10, 14, 16]. Accordingly, it is necessary to 
improve both experimental techniques and theoretical predictions. 

Theoretical studies focused on the general mechanisms behind vacuum pair cre- 
ation in the past. As those mechanisms are quite well understood by now, a point 
is reached where it is about time to perform investigations with the aim of making 
actual predictions for upcoming experiments. One way towards a more realistic de- 
scription of the Schwinger effect is to consider more realistic electric fields. Indeed, 
there have been studies which considered the electric field of a short pulse with 
sub-cycle structure recently [34, 36, 37, 38, 39, 40]: 



Here, to denotes the laser frequency, r is a measure for the total pulse length and 
(f is the carrier-envelope absolute phase. Recently, there has been an investigation 
which included the effect of a linearly varying frequency uj{t) = ui + bt as well [39]. 




(2.24) 
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It has been shown that the momentum distribution of created particles strongly 
depends on the various parameters, most notably on the carrier phase ip. It has 
also been pointed out that these distinct signatures in the momentum distribution 
might facilitate the experimental observation of the Schwinger effect because of the 
clear discrimination between particles originating from vacuum pair creation and 
background events. Further investigations could help to optimize the various laser 
parameters for upcoming experiments. My own contributions to this research topic 
[36, 37] are summarized in Section 4.4. 

2.3.2 Enhancement of the Schwinger effect 

Another topic which raised much interest very recently is the issue of dynamical 
enhancement of the pair creation rate. As the Schwinger effect is exponentially sup- 
pressed, it might happen that the reachable intensities in upcoming high-intensity 
laser facilities might still be too small so that no clear signals from vacuum pair cre- 
ation are detected. Therefore, there have been various investigations on the issue of 
how to enhance the pair creation rate by means of dynamical effects [35, 98, 99, 100]. 

The main idea of the dynamically assisted Schwinger effect is rather simple: 
Superimposing a strong and slowly varying electric field (non-perturbative pair cre- 
ation) with a weak and rapidly varying electric field (dynamical pair creation) should 
result in a decrease of the tunneling barrier and, consequently, to a drastic enhance- 
ment of the pair creation rate [35]. Indeed, it has been shown that the combination 
of the different scales results in a subtle change in the pair creation process so that 
an enhancement of the Schwinger effect can be expected [35, 99]. 

Recent investigations confirmed this picture of dynamic enhancement by means 
of quantum kinetic methods [101, 102]. Additionally, these studies gave valuable in- 
sight in the momentum distribution of created pairs and showed that the achievable 
enhancement strongly depends on the choice of the various parameters. Accordingly, 
a focus of future research should also be laid on pulse shaping analysis in order to 
further maximize the attainable pair creation rate. 

2.3.3 Depletion of the laser field 

Most investigations of the Schwinger effect considered the electric field as classical 
background. They either totally neglected the self-induced electric field due to the 
pair creation process or treated it in mean field approximation. These investigations 
indicated that the backreaction mechanism can be well ignored in the sub-critical 
field strength regime [65, 71]. Lately, there have been first attempts to account for 
the feedback of photons in a quantum kinetic formulation as well [103]. 
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Very recently, the backreaction issue beyond mean field approximation has at- 
tracted much interest [104, 105, 106, 107]. These investigations indicate that already 
a single electron-positron pair might lead to an avalanche-like electromagnetic cas- 
cade restricting the attainable laser intensities in upcoming high-intensity laser ex- 
periments: The quick acceleration in the laser field could lead to the emission of hard 
photons which in turn create electron-positron pairs via dynamical pair creation as 
already observed in the SLAC-E144 pair production experiment [108]. Accordingly, 
the evolution of a QED cascade might result in a total destruction of the laser field 
even in the sub-critical field strength regime in the last resort. It is surely worth to 
further investigate this issue because of its great importance for future high-intensity 
laser experiments. 

2.3.4 Schwinger effect in space- and time-dependent electric fields 

Investigations of the Schwinger effect mainly focused on the temporal dependence of 
the electric field E(t) in the past. It is one main result of these studies that short time 
scales yield an enhancement of the number of created particles [19, 70, 109[. On the 
other hand, there have been investigations of the Schwinger effect in spatially varying 
electric fields E(x) as well. These studies indicated that spatial inhomogeneities tend 
to decrease the pair creation rate as distinguished from temporal inhomogeneities 



This diminishment due to spatial variations can be easily understood: Vacuum 
fluctuations have to delocalize in order to gain energy and turn into real particles, 
with the Compton wavelength Ac being the characteristic delocalization scale. As 
soon as the work done by the electric field over its spatial extent is smaller than 
twice the electron mass: 



the vacuum fiuctuations cannot acquire enough energy from the electric field and 
the pair creation process terminates. 

Accordingly, there might be two opposing effects in a space- and time- dependent 
electric field E(x, t): Enhancement of the pair creation rate by temporal variations 
on the one hand and decrease of the pair creation rate by spatial variations on the 
other hand. In fact, to the best of my knowledge there have not been any first 
principle investigations on the infiuence of space- and time-dependent electric fields 
on the Schwinger effect so far. Hence, the overall effect of the interplay between the 
different scales is by no means clear. 



[22, 23, 109]. 




(2.25) 
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2.4 Objectives 

I tried to highlight some current research topics which seem of utmost importance in 
the previous section. Some of these issues wiU in fact be touched upon in this thesis, 
which is dedicated to further investigate the Schwinger effect in inhomogeneous 
electric fields. In the first part, I focus on the Schwinger effect in time-dependent 
electric fields E(t): 

• I present the close similarity between the equal-time Wigner formalism and the 
quantum Vlasov equation for the case of a time-dependent electric field E(t). 
Accordingly, it has to be emphasized that the equal-time Wigner formalism is 
the natural generalization of the quantum Vlasov equation in the presence of 
space- and time-dependent electric fields E(x, t). 

• I derive an analytic expression for the equal-time Wigner function for the 
exactly solvable cases of a static as well as a pulsed electric field. Additionally, 
I discuss some features of the Schwinger effect in a pulsed electric field with 
sub-cycle structure. It has to be stressed that this investigation represents an 
important step towards a more realistic description of vacuum pair creation 
in upcoming high-intensity laser experiments. 

As already mentioned, little is known about the Schwinger effect in space- and time- 
dependent electric fields E(x, t). Therefore, it is about time to shed some light on 
this issue. Accordingly, the following topics are covered in the second part: 

• I perform an ab initio simulation of the Schwinger effect in a simple space- 
and time-dependent electric field E(x, t) in order to better understand the 
interplay between spatial and temporal variations. Most notably, I will show 
in which way observable quantities such as the number of created particles or 
the momentum distribution change. For computational reasons, this study is 
performed on the basis of QED in 1 + 1 dimensions. 

• I adopt various approximations (local density approximation, derivative ex- 
pansion) and compare their outcome with the results of the full numerical 
simulation. These investigation will give an estimate on the range of validity 
of the different approximations. 



Chapter 3 

Schwinger effect in phase space 



In this chapter I give a brief review on the Wigner formalism upon which the sub- 
sequent investigations are based. In Section 3.1, I focus on the covariant Wigner 
formalism and derive the equation of motion for the covariant Wigner function 
[81, 83, 85]. In Section 3.2, I present the connection to the equal-time Wigner for- 
malism and discuss some of its properties in more detail [31, 33, 85, 86, 87, 92, 93]. 

3.1 Covariant Wigner formalism 

The quantum field theory of electromagnetisni is based on the QED Lagrangian:^ 

Ci^,^,A) = ^[^-f^'[i^^-eA^]^-{[i^^ + eA^]^)Jf'^]-m^^-\F^'''F^, . (3.1) 

The equations of motion for the Dirac field ^'(x), the adjoint field ^'(x) as well as 
the photon field A'^(x) are given by: 

7^[ia^-eA^(x)]*(x) -m*(x) = 0, (3.2) 

[id^ + eA^ix)] *(x)7'^ + m*(x) = , (3.3) 

and 

d,F>^'^{x) = imx),Y^{>,)] . (3.4) 

Note that these equations have to be understood as operator equations in the Heisen- 
berg picture at this point. 

Due to the fact that I focus on pair creation in an external electromagnetic field, 
the following simplification is made: The field strength tensor F^^'^{x) is not treated 
as an operator but rather as a c-number. It is clear that this mean field approxima- 
tion is only justified in the strong-field regime as it suffers from several shortcomings 
such as the neglect of radiative corrections or final state interactions. On the other 
hand, this mean field treatment results in the truncation of the BBGKY hierarchy 
for the Wigner function at the one-body level without further approximations. 

^Note, that I use the hermitian Dirac Lagrangian which treats the oppositely charged Dirac 
fields symmetrically. As usual, the field strength tensor is defined as F^'^{x) = d'^A'^{x) — d'^ A^{x). 
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3.1.1 Covariant Wigner operator 

The commutator of two Dirac fields ^ [^'^(xi), ^'q(x2)] is tlie starting point for the 
following considerations. As this density matrix is not gauge invariant under local 
U{1) gauge transformations: 

4(x) = A,{x) - a^A(x) *'(x) = vI/(x)e-^W , (3.5) 

one has to introduce a Wilson-line factor U{A;x,y): 



U{/K; X, y) = exp 



/ x+y/2 

ie I dz ■ A(z) | , (3.6) 

V 



x+y/2 \ 
e j dz A(z) , 
x-y/2 / 



in order to allow for a gauge invariant definition of the covariant Wigner operator'? 

WS(X,P) dVP-^^^(A;x,y) [^^(x- X),^^(x+ X)] , (3.7) 

with X = 5(xi +X2) being the center of mass coordinate and y = X2 — xi the relative 
coordinate.^ As already mentioned, the necessity of a gauge invariant definition is 
rather important for an interpretational reason: The Wigner function should coin- 
cide with the corresponding distribution function in its classical limit. Accordingly, 
it should be an observable quantity and as such it must not be gauge dependent. 

Note that the path-ordering in U (A; x, y) has been dropped in comparison to 
Eq. (2.20) as yl^(x) is treated as a classical mean field now. Additionally, it has 
to be stressed once more that the integration path from x — y/2 to x + y/2 in the 
Wilson-line factor is not unique. However, requiring a proper physical interpretation 
of p as kinetic momentum forces the integration path to be chosen along the straight 
line between the end points. Accordingly, VV''^)(x, p) can be written as: 



h;(4)(x,p) = ly ^^y^iP-yei-!%^^^i-+^yyy [^(x_ Z),vi,(x+ x)] . 



(3.8) 



It can be seen explicitly that the covariant Wigner operator transforms like an 
ordinary Dirac gamma matrix under Hermitian conjugation: 

[H;(4) (x, p)]t = /vv(4) (x, p)70 . (3.9) 



^Note that this definition differs from Eq. (2.19): Taking the commutator i ['J'^(xi), 'I'a(x2)] 
instead of $/3(xi)4'q(x2), the symmetrized electromagnetic current will be obtained later on. 
^Spinor indices will be suppressed in the following to simplify notation. 
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3.1.2 Equation of motion 

The equation of motion for the covariant Wigner operator can be calculated by 
taking advantage of the Dirac equation Eq. (3.2) - (3.3). Similar to the Dirac 
equation, one obtains in fact two equations which are the adjoints of each other:'^ 

[iD^(x,p)-in^(x,p)]7^W(^)(x,p) =-imW^^\x,p), (3.10) 
[iD^(x,p)+zn^(x,p)]W(4)(^^p^^M, ^ imW(4)(x,p), (3.11) 

with -D^(x, p) and n^(x, p) denoting the following self-adjoint operators:^ 

1 

77^(x,p) = cf;-e j\_AiF^,{x- 1^8^)8;, (3.12) 

1 

n^(x, p) =p^- ie /'^ dee^;.(x - iid^)d'; . (3.13) 

In order to switch from the covariant Wigner operator VV*^'^)(x, p) to the covariant 
Wigner function W^^-*(x, p), one takes the vacuum expectation value:^ 

W(4)(x,p) = {fl\W^^\x,p)\n) . (3.14) 

The corresponding equation of motion for the covariant Wigner function is accord- 
ingly found by taking the vacuum expectation value of Eq. (3.10) - (3.11). 

Note that the mean field approximation, treating F^'^{x) as a c-number instead 
of an operator, becomes apparent at this point: If the field strength tensor was still 
considered as an operator, the equation of motion for the covariant Wigner function 
{n\W^'^^n) would couple to a two-body term {n\Ft"'W^^'^\n) . As the equation of 
motion for the two-body term depends on a three-body term and so forth, this would 
generate the infinite BBGKY in the end. However, as the field strength tensor is 
in fact taken to be a c-number, the BBGKY hierarchy is already truncated at the 
one-body level. 

It seems advantageous to either add or subtract Eq. (3.10) from Eq. (3.11) for 
later use. Upon taking the vacuum expectation value, this yields two new self-adjoint 



The detailed calculation can be found in Appendix A.l. 
*The differential operators with respect to x and p are defined so that: 

d'!'. = and dt^ = as well as 9f, = and dt^ = . 

^As the field operators are considered as Heisenberg operators throughout, \Q) has to be un- 
derstood as Heisenberg vacuum as well. 
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equations of motion for the covariant Wigner function: 

lDf,[Y' -in^{7'',W(^)} = -2imW^^^ , (3.15) 

iD^{7^>V(4)}-i^^ [7^ = . (3.16) 

Along with the equation of motion for the mean electromagnetic field which results 
from Eq. (3.4), a self-consistent differential equation system is finally obtained: 

d^F^^'^ix) = (0|§ [*(x),7'^*(x)] \n) = ej-^ tr[7'^W('^)(x,p)] . (3.17) 

3.1.3 Spinor decomposition 

As the covariant Wigner function transforms as a Dirac gamma matrix, one may 
decompose W'-'^-' (x, p) in terms of the Dirac bilinears. As a matter of fact, this defines 
16 covariant Wigner components W(x, p) with specific transformation behavior:^ 

W(4)(x, p) = i [S + i75F + t'^V^ + 7^75 + a^'^T^,] . (3.18) 

Note that all W(x, p) can be chosen to be real because of Eq. (3.9). Under or- 
thochronous Lorentz transformations they transform as scalar S(x, p), pseudoscalar 
P(x, p), vector V^(x, p), axialvector A^(x, p) and tensor T^[/(x, p), respectively. 

In order to decompose Eq. (3.15) - (3.16) in the basis of Dirac bilinears, one 
has to calculate both the commutator and the anticommutator between 7^ and the 
Dirac bilinears. This gives:^ 





{7^-} 


[7^■] 


I 


27/^ 





75 





27/^75 




2gt^- 








2<'^75 






2i{g'"'Y - gPPY) 



(3.19) 



'^Alternatively, it is also possible to first multiply Eq. (3.10) from the left and the right side by 
a factor 7°. In this case, the resulting equation of motion becomes symmetric under the exchange 
of commutators [■,•]_ with anticommutators [•, -J^: 

[7V, ^ - tii, [7%^ w^^h°] ^ = -i^ [7°, yv^^h"] ^ ■ 

^The components of the covariant Wigner function will be denoted by upper-case fonts W(x, p), 
whereas the components of the equal-time Wigner function will be denoted by lower-case fonts 
w(x, p, t) in the following. 

^Here, the usual convention e"^^-^ = 1 for the totally antisymmetric tensor is used. 
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Accordingly, the terms appearing in Eq. (3.15) - (3.16) read: 

[Y, W(^) } = i [t^S + W + y'^^a^pK, - e^'-P^^^^^T.p] , (3.20) 
[jt^, H;(^)] = i [i7^75P - icT^'V, + 75A^ + 2igf'''YT.p] ■ (3.21) 

The equations of motion for the covariant Wigner components are then found by 
equating the coefficients of the individual Dirac bilinears. On the one hand, this 
results in a set of inhomogeneous equations originating from Eq. (3.15): 

WVf,, = mS , (3.22) 

D^A^, = 2mP , (3.23) 

n^S - Dj.T"'' = mV^ , (3.24) 

D^P - 2e^"'P''U„Tp^ =-2mA^ , (3.25) 

(L^'V - D^V) + 2e>"'P''UpA^ = AmTf"" . (3.26) 



On the other hand, a set of homogeneous equations is derived from Eq. (3.16) 
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Df'Vp = 0, (3.27) 

n^A^ = 0, (3.28) 

DfS + 4U^T''>' = 0, (3.29) 

n^^P + igM^p^IP^Tp^ = 0, (3.30) 

- U-'V) - ^ef^P^DpA^ = . (3.31) 



3.2 Equal-time Wigner formalism 

The covariant Wigner formalism provides a clear theoretical description, however, 
the experimental situation rather suggests a non-covariant formulation in terms of 
an initial value problem: Accordingly, covariance should be given up in favor of a 
clearer and simpler interpretation. In fact, there are two different but equivalent 
ways to obtain an equal-time Wigner formalism: One can either derive the equation 
of motion for the equal-time Wigner function W*-^-' (x, p, t) in a non-covariant way 
from the very beginning or start with the covariant formulation and perform an 
energy average afterwards. ""^^ 

^"The set of equations Eq. (3.22) - (3.31) is sometimes called VGE equations in the literature, 
named after D. Vasak, M. Gyulassy and H. T. Elze [83]. 

^^The equal-time Wigner formalism has first been derived in a non-covariant way from the very 
beginning [31] whereas the connection between the covariant and the non-covariant approach has 
been shown only later ]85, 86, 87]. Note that the second approach is adopted in the following. 
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3.2.1 Energy average 

The equal-time Wigner function is defined as the energy average of the covariant 
Wigner function: 

W(3)(x,p,t)^ I ^wW(x,p). (3.32) 
Plugging in the explicit expression for the covariant Wigner function, this reads: 

W^^nx,?,^ = IJ dV"'P-^e-'^^-V^''^'^^"^^'''*^''^(f^l[^(x- |,0,*(x+ 

(3.33) 

with t = . The equal-time Wigner function can then be decomposed into its Dirac 
bilinears in complete analogy to Eq. (3.18):^^ 

>V<^^ (x, p, t) = 1 [s; + i75p + + 7^753^ + a^^Hf,^] . (3.34) 

Note that the equations of motion for the covariant Wigner components cannot 
be set up as an initial value problem by fixing an asymptotic field configuration at 
— > —00. This is in contrast to the equal-time approach which will indeed allow 
for a formulation in terms of an initial value problem. -"^^ It has to be emphasized, 
however, that dynamical information on W(^) (x, p) is lost by taking the energy aver- 
age. As a consequence, in order to recover all the information which is thrown away 
by averaging out the off-shell behavior of VV'-'*^(x, p), one would have to consider not 
only the equations of motion for w(x, p, t) but also for all higher energy moments: 

wW(x,p,t) = J ^p^W(x,p) . (3.35) 

Accordingly, the equal-time Wigner formalism would have to involve the equations 
of motion for all energy moments w["'l(x, p,i) in order to provide a complete de- 
scription of the system. This would, however, result in an infinite hierarchy of equa- 
tions again. This hierarchy is in fact truncated at the level of the equal-time Wigner 
components w(x, p, t) for an electromagnetic field which is treated as classical mean 
field. Accordingly, the equal-time Wigner formalism in mean field approximation is 
complete in the sense that knowledge of w(x, p, t) suffices to recover all information 
contained in the covariant Wigner function W^'^\x,p) [92, 93]. 

^^The equal-time Wigner components are defined by the energy average as well: 

w(x,p,i) = /^W(x, p). 

^^The simplest choice of initial condition, which will be used later in this thesis as well, is to 
take the vacuum value of the equal-time Wigner components Wvac(x, p, tvac), with t 

vac ^ OG . 
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3.2.2 Equation of motion 



The derivation of tlie equation of motion witliin tlie equal-time Wigner formalism is 
based on the equations of motion for the covariant Wigner components. Taking the 
energy average of Eq. (3.22) - (3.31), it turns out that the subset of equations which 
does not contain any explicit dependence on po via Uq (x, p) trivially translates into 
time evolution equations^^ for the equal-time Wigner components:^^ 



Dtp 
Dt vo 
Df ao 



- 2n- ti 
+ 2n ■ t2 



D ■ V 
D a 

D vo + 2n X a 

DfSL + D ao + 2n X w 

Dt ti +1) xt2 + 2Usi 

A*2 - D X ti - 2np 





-2m ao 


2m p 
-2m ti 


2m V 




(3.36) 
(3.37) 
(3.38) 
(3.39) 
(3.40) 
(3.41) 
(3.42) 
(3.43) 



with the vectors ti(x, p, t) and t2(x, p, t) being defined in terms of the antisymmet- 



ric tensor t 



.16 



tl 



2t*ei 



and 



t2 



(3.44) 



The pseudo-differential operators Df {x,p,t), D(x, p, t) and Tl{x,p,t), which are 
obtained by taking the energy average of Eq. (3.12) - (3.13), are given by: 



A(x,p,t) = I 
D(x,p,t) = Vx 
n(x,p,i) = p- 



^deE(x + ieVp,i) ■ Vp , 

2 
L 

[ d^B{x + i^Vp,t) X Vp , 

2 
L 

\ daB(x + ieVp,i) X Vp 



(3.45) 
(3.46) 
(3.47) 



As a matter of fact, the equations of motion for the covariant Wigner components 
Eq. (3.22) - (3.31) give rise to still another set of constraint equations, resulting 



^■^The set of equations Eq. (3.36) - (3.43) is also called BGR equations in the literature, named 
after I. Bialynicki-Birula, P. Gornicki and J. Rafelski [31]. 

^^One must assume that the electromagnetic field can be expanded in a Taylor series with respect 
to the temporal coordinate and that the covariant Wigner components and all its derivatives with 
respect to po vanish for |po| — > oo. Details can be found in Appendix A. 2. 

'^^This is in analogy to the definition of the electric and magnetic field in terms of _F'"'(x). 
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from the energy average of the subset of equations which does show an exphcit 
dependence on pq via IIo (x, p) : 







- iD ■ tl 




m vq , 


(3.48) 
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+ 5D ■ t2 
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(3.49) 
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(3.50) 
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-n vo = 


, 


(3.52) 


all] 


+ Hia 


+ 5D X V 


— n ao = 


— m ±2 , 


(3.53) 




+ Hiti 


+ iDs 


-n X t2 = 


, 


(3.54) 




+ Hits 


- iDp 


+n X ti = 


— m Si . 


(3.55) 



This system of equations describes the relation between the first energy moments 
w[il(x, p, t) and the equal-time Wigner components w(x, p,t). Moreover, the oper- 
ator Ilt{x,p,t) is given by: 

1 

nt{x,p,t)=ie dCeE(x + ieVp,t)-Vp . (3.56) 

The separation into transport and constraint equations continues for all higher 
energy moments as well: Multiplying the equations of motion for the covariant 
Wigner components Eq. (3.22) - (3.31) by a factor and subsequently taking the 
energy average results in a set of transport equations for w["l(x,p,i) and a set of 
constraint equations. As mentioned in the previous section, this finally gives an 
infinite hierarchy of equations for all energy moments. 

In mean field approximation, however, this infinite hierarchy is truncated at the 
level of the equal-time Wigner components w(x, p, t). This means that the transport 
equations for w["l(x, p, t) are identically fulfilled once w(x, p, t) obey the equations 
of motion Eq. (3.36) - (3.43).^^ Accordingly, in order to calculate the higher energy 
moments w^"] (x, p, t) one does not have to solve the infinite hierarchy but it suffices 
to deduce them from w(x, p, t) via the constraint equations. 

Due to the fact that the equal-time Wigner components w(x,p,i) encode the 
information of all higher energy moments w["l(x,p,i), the equal-time Wigner for- 
malism based on the equations of motion Eq. (3.36) - (3.43) is called complete. 

^^For the first energy moments w'^' (x, p, t) this is demonstrated in Appendix A. 3. The general 
proof can be found in [93]. 
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3.2.3 Observables and conservation laws 

It is one objective of this tliesis to calculate observable quantities such as the mo- 
mentum distribution or the total number of created particles for the Schwinger effect 
in space- and time-dependent electric fields E(x,f) later on. 

As a matter of fact, the observables in this theory are defined via Noether's 
theorem and therefore correspond to conserved quantities, most notably the total 
charge Q, the total energy S, the total linear momentum P, the total angular 
momentum M as well as the Lorentz boost operator K. 

3.2.3.1 Electric charge 

The invariance of the QED action under global U{1) transformations: 

$'(x) = *(x)e*" , (3.57) 
results in the electromagnetic current: 

j^(x) = |[*(x),7^^(x)] , (3.58) 

with: 



d^J^ix) = . (3.59) 

The vacuum expectation value of this electromagnetic current can in fact be ex- 
pressed in terms of the covariant Wigner function W^^\x, p): 

{n\f{x)\n) = e| ^ tr [7'^W(4)(x, p)] , (3.60) 

so that the space integral of its 0-component, which corresponds to the total electric 
charge Q, is constant in time. Expressing it in terms of the equal-time Wigner 
components, it is given by: 

Q = J d^x{n\f{x)\n) = e J drvo(x,p,t) , (3.61) 

with dr = d'^a:d'^p/(27r)'^ being the phase-space volume element. 
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3.2.3.2 Energy and linear momentum 

Invariance of the QED action under space-time translations results in the canonical 
energy-momentum tensor T'"^(x), with:^* 

a;^r^^(x) = . (3.62) 

The energy-momentum tensor T'*'^(x) is composed of two parts. The electromagnetic 
part, which is in fact symmetric and gauge-invariant, reads: 

TZi^) = \9^''F^^{x)F^p{x) - g''"F^/^(x)F„^(x) , (3.63) 

whereas the gauge-invariant fermionic part, which treats the oppositely charged 
Dirac fields symmetrically, is given by [112]:^^ 

T^'{x) = \ [^'(x),7'^K - eA''(x)]*(x)] - \ {W^ + eA''(x)]f'(x),7"*(x)] . (3.64) 

It is again possible to express the vacuum expectation value of the fermionic part 
in terms of the covariant Wigner function W''^-'(x, p): 

(0|Tr(x)|0> = I [7^W(4Hx,p)] . (3.65) 

As the energy-momentum tensor is conserved, the space integral of the 0-coniponents 
are again constant in time. In fact, this defines the total energy £ and the total 
linear momentum P: 



d3x(0|T°°(x)|0) and P = / d^x (0|r°Xx)e,|Q) . (3.66) 



^*The canonical energy-momentum tensor T^" {y) obtained from Noether's theorem is neither 
symmetric nor gauge invariant. It is, however, always possible to find a remedy by adding a total 
derivative term which is antisymmetric under the exchange of and A [110, 111]: 

r^''(x) = T""^(x) + ajs^''^(x) = r^''(x) + aj[F"^(x)/i''(x)] . 

^^The fermionic part of the energy-momentum tensor [y) is still not symmetric. A symmetric 
definition would be given by the linear combination: 

r£:,,„(x) = i[r/,"'(x) + r^''(x)] . 

Note, however, that a symmetric definition is not absolutely necessary for the purpose of defining 
the total energy E and the total linear momentum P. 
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Expressing S in terms of the equal-time Wigner components, this gives:^" 
£ = j dr[mffi(x,p,i) +p- v(x,p,i)] + i yd^a; [|E(x,i)|2 + |B(x,i)n , (3.67) 



whereas P reads: 



j drpvo(x,p,t) + j d^xE(x,t) X B(x,t) . 



(3.68) 



3.2.3.3 Angular momentum and Lorentz boost 



Invariance of the QED action under homogeneous Lorentz transformations results 
in the generalized angular momentum tensor J^'^'*'(x), with: 

dy^^ix) = . (3.69) 

Treating the oppositely charged Dirac fields symmetrically, it is given by [110]:^"'^ 

J'-^(x) = x''Tt/^{x) - x''T,7„(x) + i [*(x), {7^ ^} *(x)] . (3.70) 

The first two terms describe the orbital part whereas the last term originates from 
the spin of the Dirac fields. The vacuum expectation value of the orbital part can 
be taken from the previous section, whereas the spin contribution is given by:^^ 

= 2 / (^t^- [{7^^}wW(x,p)] . (3.71) 



^"in order to calculate the fermionic part of the total energy, the expression for the first energy 
moment Wq^'(x, p,t) in Eq. (3.50) is needed: 



dx {n\TS {x)\n) = J dr^j,"' (x, p, t) . 

As an integration over the whole phase space is performed, this can be simplified to give: 
J d^x (fi|rD°(x)|n> = J dT[ms + n • w - Htwo] = J dT[ms + p ■ w] . 

^^Note, that the symmetric energy-momentum tensor Ts'^m(x) has to be used in this expression. 
As one is forced to use the symmetric energy- momentum tensor, the fermionic part is given 

by: 

mZyn-Mm = 5 / (i^ tr [{p^j^ +p^Y)w^'\^, P)] • 

Even though this expression differs from Eq. (3.65), it turns out that the integration over the whole 
space gives the same result: 

|d^:r(fi|r«''(x)|n> = I d'^x{n\T^%,,,{x)\Q) , 

such that the observable quantities Eq. (3.67) and Eq. (3.68) remain unchanged. 
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As the generalized angular momentum tensor is conserved, the space integral of the 
0-coniponents are again constant in time. In fact, this defines the total angular 
momentum M and the Lorentz boost operator K: 

M = J d^x{n\^eijkJ°^^{^)ei\n) and K = y d^x (0| J°°'(x)ei|0) . (3.72) 

Expressing these quantities in terms of the equal-time Wigner components, M is 
given by: 

M= y dr [x X pvo(x,p,i) - ia(x,p,i)] d^.xx X E(x,t) X B(x,t) , (3.73) 
whereas K reads: 

K = tP- J drx[ms(x,p,i) +p- v(x,p,t)] - ^J d3xx[|E(x,i)|2 + |B(x,t)|2] . 

(3.74) 



Chapter 4 

Schwinger effect for E(t) 



This chapter is dedicated to the quantum kinetic formulation of the Schwinger ef- 
fect in spatiaUy homogeneous, time-dependent electric fields E(t). In Section 4.1, 
I briefly review the flrst principle derivation of the quantum Vlasov equation in its 
nowadays widely used form [29, 30, 65]. In Section 4.2, I show that the equal-time 
Wigner formalism for spatially homogeneous, time-dependent electric fields is in fact 
totally equivalent to the quantum Vlasov equation. In the subsequent Section 4.3, 
I derive analytic expressions for both the one-particle distribution function and the 
equal-time Wigner components for the exactly solvable cases of a static as well as a 
pulsed electric field [33]. In Section 4.4, I eventually discuss the Schwinger effect in 
a pulsed electric field with sub-cycle structure on the basis of a numerical simulation 
of the quantum Vlasov equation [36, 37[. 



4.1 Quantum Vlasov equation 

As already mentioned, it was only in the mid 1990s that a first principle derivation 
of a quantum kinetic equation accounting for electron-positron pair creation in the 
presence of a spatially homogeneous, time-dependent electric field E(i) has been 
given. Since then a huge number of studies for various types of time dependencies 
has been performed and achieved in fact significant progress [36, 37, 57, 66, 67, 68, 
69, 70, 102, 103]. 

This formalism is based upon the QED Lagrangian Eq. (3.1), with the field 
strength tensor treated as as c- number again. ^ Accordingly, the equation of motion 

^In this mean field approximation, it is assumed that the magnetic field vanishes whereas the 
electric field is time dependent and points along a given direction: 

B{t) = and E(t) = E{t)e3 . 

It is convenient to represent this type of electromagnetic field in temporal gauge Ao{x) = 0. 
Accordingly, the electric field is derived from the vector potential: 

E(t) = -iA(t) with A{t) = A{t)e3 . 
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for the Dirac field Eq. (3.2) reads: 

[Vx-ieA(i)] -m] *(x,i) = . (4.1) 

A decomposition of tire Dirac field into its Fourier modes is indicated as the vector 
potential is spatially homogeneous: 

(2^)3 



*(x,i)= / 7;^e^''-"^(q,t) , (4.2) 



so that the corresponding equation of motion for the Fourier modes reads:^ 

[i/|_^.^(q,t)_^]^(q,i)=0. (4.3) 

This equation is solved by introducing the following ansatz: 

V'(q, t) = [i7° I - 7 ■ 7r(q, i) + m] 0(q, t) , (4.4) 

so that the first order differential equation for t) is transformed into a second 
order differential equation for (f){q,t):^ 

'g, + LO^q, t) + zeii;(t)7°7'] ^(q, t) = . (4.5) 

It is convenient to expand 0(q, t) in a basis^ spanned by the eigenvectors of 7'^7'^: 

4 

ct>{q,t)=J29r{ci,t)Rr . (4.6) 

r=l 



^The following notation is used throughout: The variable q is a canonical momentum, whereas 
7r(q, t) = q — eA{t) denotes the time-dependent kinetic momentum on a trajectory. Moreover, this 
has to be clearly distinguished from the kinetic momentum in phase space p. 

^It is convenient to introduce the following notation for later use: The transverse energy squared 
is defined as = + Q± so that the total energy squared is given by: 

aj^(q, t) = + 7r^(q, t) = + ^3(93,*) ■ 

''The Dirac bilinear 7''7'^ = — icr"^ is proportional to the Lorentz boost operator along the 
electric field direction K3 = 50" ■ The four corresponding eigenvectors R,-: 

K3Rr={l,2} = +5-R7={l,2} and -K"3-Rr-={3,4} = -5-Rr={3,4} 

form a complete basis. As the commutator [_R'3,A/3] = vanishes, where AI3 = icr^^ denotes 
the spin operator along the electric field direction, one may further characterize the eigenvectors 
according to their intrinsic spin: 

M3_Rr={l,3} = +5-Rr={l,3} and M3_R,.= {2,4} = — 5-Rr={2,4} ■ 
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Inserting this ansatz into Eq. (4.3), it turns out that each gr{q,t) obeys the equation 
of motion of a time-dependent harmonic oscihator with an additional complex term: 



^ + w^(q,^) 



ieE{t) gr{q,t) 



92 



+ uj^{q,t) -ieE{t) gr{q,t) 



with r = {1,2} , (4.7) 
with r = {3,4} , (4.8) 



which are in general not exactly solvable; exceptions are the constant electric field 
E{t) = Eq and the pulsed electric field E{t) = EQsech'^{^) to be discussed in more 
detail later on. Note that Eq. (4.7) - (4.8) are second order ODEs, each of them pos- 
sessing two linearly independent solutions gi^\q,t), so that (j){q,t) is composed of 
eight independent contributions in total. ^ In contrast, the Dirac equation Eq. (4.3) 
allows merely for four linearly independent solutions. This redundancy is removed 
by choosing only one complete set of solutions, either the one with r = {1,2} or 
the one with r = {3,4} [28]. Choosing the first set, for instance, Eq. (4.4) can be 
written as a linear combination: 



z 

^P{q,t)=Y^ [ur{q,t)ai+\q)+Vr{-q,t)al-\-q) 



r=l 



(4.9) 



with ai^\q) being the coefficients of the corresponding terms. Moreover, the spinors 
are given by:^ 



,,(q, t) = [i7°| - 7 ■ 7r(q, t) + m] g^+^ (q, t)Rr , 
,(-q, t) = [^7°! - 7 ■ 7r(q, t) + m] g^-^ (q, t)Rr . 



(4.10) 
(4.11) 



The coefficients are interpreted as anticommuting creation/annihilation operators 
for particles and antiparticles upon canonical quantization: 



q) — > ttriq) and 



f'J(-q) , 



(4.12) 



with the only non- vanishing anticommutators being given by: 

{a,(q),at(q')} = {6,(q) , 6t(q')} = {2nf 5rs5{q- c{) . 



(4.13) 



^In the limit of a vanishing electric field, they turn into the positive and negative energy plane 
wave solutions: 

3r|^L(q, t) ~ e'''""*'''' with cj^(q) = + . 

®The linearly independent solutions gi*'(q, and gi^\ci,t) are identical due to the uniqueness 
theorem for ODEs: 

gJ±'(q,t)=gfHq,t)=gf±'(q,t) . 
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It has to be emphasized that a particle interpretation of the field quanta is only 
possible in the case of plane wave solutions.^ The actual mode functions g^^\q^,t), 
however, are no plane waves once an electric field is present so that a particle 
interpretation is not feasible. A further consequence of a non-vanishing electric field 
is that the Hamiltonian operator achieves off-diagonal elements, accounting for the 
possibility of pair creation as well as pair annihilation. 

The off-diagonal Hamiltonian operator can be diagonalized by performing a uni- 
tary non-equivalent change of basis to a quasi-particle representation via a time- 
dependent Bogoliubov transformation: 

Ar{q,t) = a(q,t)a,(q)-/?*(q,t)4(-q) , (4.14) 
Bl{-q,t) = (3{q,t)ar{q)+a*{q,t)bl{-q) . (4.15) 

In order to be a canonical transformation, the Bogoliubov coefficients a(q, t) and 
/3(q, t) must fulfill: 

|a(q,i)l' + l/?(q,i)P = l • (4-16) 

The creation/annihilation operators, even though they are now considered as being 
time-dependent, still fulfill the equal-time anticommutation relations: 

{A,(q,t),4(q',0} = {Br{q,t), Bl{q',t)} = (27r)3<5,,5(q - q') . (4.17) 

In the vacuum, when there is no electric field present, the Bogoliubov transformation 
is trivial and the two operator bases coincide: 

a(q,ivac) = l and /3(q,tvac) = with tvac ^ -oo . (4.18) 

In fact, the time-dependent Bogoliubov transformation is equivalent to expanding 
the Fourier modes 'i/;(q, i) in a different basis: 

2 

i:{q,t)=Y,[Ur{q,t)Ar{q,t) + Vr{-ci,t)Bl{-q,t)\ . (4.19) 

r=l 

As a matter of fact, any complete basis of spinors Ur{q,t) and Vr{q,t) would do 
this job. The adiabatic basis, however, which is chosen in close analogy to the basis 



''The spinor normalization is chosen according to uj.(q, t)us(q, t) = vl{q,t)vs(q,t) = Sra, so 
that the properly normahzed positive and negative energy plane wave solutions are given by: 

ffvac(q,t) = 
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of spinors Wr(q, i) and Ur(q, i), respectively, fits especially well: 

Ur{ci,t) =[j°Uj{(l,t)--f7T{(l,t)+m]G^+\ci,t)Rr, (4.20) 
K(-q,i) = [-7°c^(q,t)-7-7r(q,t)+m]G(-)(q,i)i?, , (4.21) 

with the adiabatic mode functions G'^^^(q, i) being defined as:^ 

.,. gTie(q,to,t) ft 

G(^)(q,0= , ^ .^ . .m ^^^^ e{q,to,t) = u;{q,t')dt' . 

V2a;(q,t)[w(q,t) =F7r3(g3,t)J ^to 

(4.22) 

The lower bound to of the dynamical phase G(q, t;fo) is not fully determined as it 
only fixes an arbitrary phase. Note that the adiabatic mode functions G^^\q,t) 
behave like plane waves in the limit of a vanishing electric field. ^ 

The one-particle distribution function J-{q, t) is defined as the instantaneous 
quasi-particle number density for a given canonical momentum q. As quasi-particles 
with different values of r behave identically when there is no magnetic field present, it 
is in fact convenient to include the sum over both values r = {1, 2} in the definition:^'' 

^(,, , H,n ± '"l^'fa'y-'"'"!"' ^ OP . (4.23) 

V—>oo '—^ V 
r=l 

with vacuum initial conditions for iyac —oo being assumed: 

(0|4(q)a,(q)|0) = (0|6j(q)6,(q)|0) = . (4.24) 

Accordingly, the knowledge of /3(q, t) allows for the calculation of t). Note that 
the difi'erent representations Eq. (4.10) - (4.11) and Eq. (4.20) - (4.21), respectively, 
translate into a relation between the Bogoliubov coefficients: 

a(q, t) = ie^G^-^ (q, t) [| - ic^(q, t)] g^+^ (q, t) , (4.25) 
/3(q, t) = -ze^G^+^ (q, i) [| + t)] (q, t) , (4.26) 



*In fact, this choice is very convenient as I7r(q, t) and Vl (q, t) coincide with it7 (q, t) and fr(q, t), 
respectively, in the hmit of a vanishing electric field. 

®This is the deeper reason why the one-particle distribution function J^(q, t) can be interpreted 
as the momentum distribution of real particles at asymptotic times when the electric field vanishes. 

^"Here, V denotes the configuration space volume, which has to be taken to infinity. This factor 
is necessary to cancel the divergence which appears in the anticommutators: 

{a,.(q),at.(q)} = {&,.(q) , foj, (q)} = (2^)-^5(0) = V , 

Alternatively, one could derive the quantum Vlasov equation in a finite box from the very beginning 
so that q can take only discrete values. In this case, the limit V ^ oo would be taken at the very 
end of the derivation. 
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so that their time derivatives form an ODE system: 



«(q,i)= iQ(q,i)/3(q,0e''®('''*°'*^ , (4-27) 
/3(q, t) = -iQ(q, t)a{q, t)e-2^e(q.*o.*) , (4.28) 

with Q(q, t) being defined as: 

Introducing the auxiliary quantity C(q, t) which describes the density of quasi- 
particle pair creation: 

V—>oo ^— ^ 1/ 
r=l 

the ODE system Eq. (4.27) - (4.28) can be rewritten in terms of the functions C(q, t) 
and J='{q, t):'^^ 

C(q, t) = - Q(q, t) [1 - ^(q, t)] e-2^e(q.*o,t) ^ (4.31) 
-^(q, t) = - Q(q, i) Re [C(q, t)e^''^(^^'^^^)] . (4.32) 

Integrating the first equation from tyac to t yields the quantum Vlasov equation for 
the one-particle distribution function .F(q, t) in its integro-differential form [30]: 



J^(q,i) = Q(q,i) J Q(q,t')[l-^(q,iO]cos[2e(q,t',t)] , 



(4.33) 



along with the vacuum initial condition .F(q, tvac) = 0. Accordingly, the pair cre- 
ation process shows non-Markovian behavior owing to memory effects and non- 
locality in time. Additionally, [1 — t)] might be interpreted as a Pauli-blocking 
factor which accounts for the correct quantum statistics. 



It turned out to be advantageous to rewrite this integro-differential equation in 
terms of a first order ODE system as the numerical solution of integro-differential 
equations is usually quite challenging [65]. To this end, two auxiliary quantities are 

^^The sum over r = {1,2} is taken again as one is only interest in the total density of quasi- 
particle pair creation. Additionally, the infinite volume limit V ^ 00 has to be taken in order to 
cancel the divergence originating from the anticommutators. 

""^^In order to derive these equations, the relation ^|/3(q, t)|^ = 2 Re[/3(q, t)/3*(q, t)] is used. 
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introduced: 



t 

Q{c\,t) = y Q(q,iO[l-^(q,0]cos[2e(q,i',t)] , (4.34) 

^vac 
t 

n{q,t) = J Q{ci,t')[l-T{q,f)]sm[2e{q,t',t)] , (4.35) 



so tliat: 



T{q,t) = Q{q,t)giq,t) , (4.36) 
g{q,t) = Qiq,t)[l-J^iq,t)]-2ujiq,t)n{q,t), (4.37) 
n{q,t) = 2Lo{q,t)g{q,t) , (4.38) 

along with vacuum initial conditions .F(q, tyac) = G{q,tvax:) = "^(Qi^vac) = 0. 

As mentioned previously, a consistent mean field description of the Schwinger effect 
requires to include the self-induced electric field due to the pair creation process: 
The charged particles, which are created by the applied electric field Eext(^) serve 
as a source of an internal electric field Eint(i)- Accordingly, one is faced with the 
problem of describing the pair creation process in the presence of the total electric 
field E(f) = Eext(0 + Eint(i) in a self-consistent way. 

The internal electric field in mean field approximation can in fact be calculated 
from the vacuum expectation value of the electromagnetic current:-"^^ 

^int(t) = -{n\f{t)\n) = -§/ d'x(n| [*(x,t),73^(x,0] |0) . (4.39) 

As a matter of fact, this expression exhibits a logarithmic UV-divergence, which has 
to be regularized and renormalized. This finally yields [65]: 



^int(t) = -2e 



^4^^(q,t) + ^Hq,t) - 
uj{q,t) ^ eE{t) ^ 4a;5(q,t) 



(4.40) 



The first two terms represent the conduction and polarization current, respectively, 
whereas the last term is introduced as a counterterm in the course of charge renor- 
malization. Together with the quantum Vlasov equation, this provides a consistent 
mean field description of the Schwinger effect in time-dependent electric fields. 



^^As I will not consider the backreaction issue in the following, I only give the final result for 
completeness' sake. Details can be found in [53, 65, 71]. 

^*Note, that Eint(t) points along the same direction as Eext(i) by virtue of cylindrical symmetry. 
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4.2 Equal-time Wigner formalism 

There are now two at first sight rather different formalisms at hand: The equation 
of motion for the equal-time Wigner function has been derived by Wigner trans- 
forming the vacuum expectation value of the fermionic density matrix whereas the 
quantum Vlasov equation results from canonical quantization of the Dirac field and 
subsequent Bogoliubov transformation to a quasi-particle representation. 

It has to be emphasized that the equal-time Wigner formalism allows for an 
appropriate description of the pair creation process in the presence of a space- and 
time-dependent electromagnetic field F^'^{x) whereas the quantum Vlasov equation 
is only valid for spatially homogeneous, time-dependent electric fields E(i) = E{t)es- 
Accordingly, the question arises whether there is a deeper connection between them 
so that the equal-time Wigner formalism can in fact be related to the quantum 
Vlasov equation. 

As already mentioned, it is an advantage of the equal-time Wigner formalism 
that it is formulated in terms of an initial value problem. The appropriate initial 
conditions for describing the Schwinger effect are provided by the vacuum value 
of the equal-time Wigner function, which can be easily calculated by explicitly 
evaluating Eq. (3.33) for the free Dirac field. The only non-vanishing equal-time 
Wigner components are then given by:^^ 

- and Vvac(p) = T-r ■ (4.41) 



a;(p) w(p) 

They are in fact spatially homogeneous and time-independent as it should be. 

Considering the equations of motion for the equal-time Wigner components 
Eq. (3.36) - (3.43) in the presence of a time-dependent electric field E(i) = E{i)&-},, 
the pseudo-differential operators Eq. (3.45) - (3.47) reduce to local ones:^'' 

A(x,p,0 I ' (4.42) 

D(x,p,i) Vx , (4.43) 
n(x,p,i) p . (4.44) 



^^The equal-time Wigner formalism in the limit of a time-dependent electric field has in fact been 
considered already at the time when the formalism was invented [31]. However, as the nowadays 
widely used form of the quantum Vlasov equation has not been derived by then [30, 65[, the deep 
connection between them has been pointed out only recently in the context of this thesis [33[. 

^^The detailed calculation can be found in Appendix A. 4. 
For a more general time-dependent electromagnetic field F'"^{t), these operators are given by: 

Dt^§l+ eE(t) ■ Vp and D ^ Vx + eB x Vp and H ^ p . 
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Moreover, all spatial derivatives vanish owing to the spatial homogeneity of both 
the electric field E{t) and the vacuum value of the equal-time Wigner components 
Wvac(p)- Accordingly, the equations of motion for the spatially homogeneous equal- 
time Wigner components w(p, t) are given by: 



As 
Dtp 
Dt vo 
Dt ao 

Dt a 
A t2 



2p ■ ti 

2p-t2 



2p X a 
2p X V 

2pfB 

2pp 





-2m ao 


2m p 
-2m ti 


2m V 




(4.45) 
(4.46) 
(4.47) 
(4.48) 
(4.49) 
(4.50) 
(4.51) 
(4.52) 



First of all, the equal-time Wigner component vo(p, t) decouples from this first order, 
partial differential equation (PDE) system. Due to the fact that vo,vac(p) = 0, this 
yields:^* 

vo(p,i)=0. (4.53) 

Moreover, the equal-time Wigner components p(p,t), ao(p, i) and t2(p,i) form a 
closed set. As their vacuum values vanish as well, one finds: 

p(p,0 = ao(p,t) =t2(p,0 =0 . (4.54) 

As a consequence, one is left with the PDE system: 

|+ei?(0^]^(p,0 = M(p)^(p,t) , (4.55) 

where ^(p, t) = {a, v, a, iti}(p, t) denotes the vector consisting of the remaining ten 
equal-time Wigner components. Moreover, M(p) is a shorthand notation for the 
10 X 10 matrix in accordance with Eq. (4.45) and Eq. (4.49) - (4.51). As Eq. (4.55) 
is in fact a hyperbolic PDE system, it can be solved by adopting the method of 
characteristics.^^ Formally, this is done by replacing the kinetic momentum in phase 
space p by the time-dependent kinetic momentum on a trajectory 7r(q, t) so that 



According to the interpretation of vo(p,it) as the charge density, this result is interpreted as 
local charge neutrality in the spatially homogeneous situation. 

^^The differential operator in Eq. (4.55), acting on any differentiable function of the independent 
variables p and t, might be viewed as a total derivative with respect to a new parameter a, assuming 
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the PDE system is transformed into an ODE system:^'^ 

^(q,t) = M(q,i)^(q,t) . (4.56) 

As ^(q, t) is a 10-component vector, it is suggested to search for an orthonormal 
basis of the corresponding 10-diniensional space in such a way that: 

10 

^(q,t) = -2J2SM, t)^M,t) . (4.57) 

Here, ei(q, t) are the basis vectors whereas £'i(q, t) denote the expansion coefficients. 
The factor —2 is chosen for later convenience. Exploiting the vacuum initial con- 
ditions, it is convenient to choose the first basis vector (Bi(q, i) in such a way that 
£^i(q, tvac) = 1 whereas all other £^i(q, tvac) = 0. As a consequence, one finds that a 
closed subset of only three basis vectors couples to the vacuum initial conditions:^^ 

M(q,f)a3i(q,i) = ei(q, t) = -Q(q, i)e2(q, i) , (4-58) 

M(q,i)e2(q,t) = 2t^(q, i)e3(q, t) e2(q,t)= Qiq,t)si{q,t) , (4.59) 

M(q,t)e3(q,t) = -2t^(q,i)<B2(q,i) e3(q,i)= 0, (4.60) 

so that the system is fully characterized by means of three expansion coefficients 
£i{q,t) with i = {1,2,3}. According to Eq. (4.56), one finds: 

£i{q,t) = -Q{q,t)£2{q,t) , (4.61) 
4(q,i) = Q{q.t)£i{q,t) - 2Lo{q,t)£s{q,t) , (4.62) 
4(q,0 = 2uj{q,t)£2{q,t) , (4.63) 

along with vacuum initial conditions £i{q, iyac) = 1 and ^^2(q) ^vac) = ^3(q) ^vac) = 0. 



that the former independent variables depend on a now: 



[±+eE{t)^]f{p,t)^£fip{a),t{a)) , 

with t{a) = t and p(a) = 7r(q, t). Accordingly, the kinetic momentum in phase space p is traded 
for the time- dependent kinetic momentum on a trajectory 7r(q, = q — eA{t), with q serving as 
an integration constant and corresponding to the canonical momentum. 

^"The following notation will be used throughout: Any function depending on p and t is under- 
stood as being defined in phase space, whereas any function depending on q and t is considered on 
a trajectroy. They are related to each other via: 

/(q-*) = /(p,t)lp^q-eA{t) and /(p, t) = /(q, i)|q-p+eA(t) • 

Note, that /(p, t) possesses only an explicit time dependence whereas /(q, t) shows both an explicit 
and an implicit time dependence via 7r(q,t). 

^^The detailed calculation can be found in Appendix A. 5. 
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Introducing the notation: 

^i(q,i) = l-^(q,t), (4.64) 

52(q,t) = g{<l,t) , (4.65) 

Ss{ci,t) = nici,t) , (4.66) 

so that J-{q, t) measures the deviation from the vacuum state, one exactly recov- 
ers the quantum Vlasov equation Eq. (4.36) - (4.38) in its diflFerential form along 
with identical initial conditions. Accordingly, one can conclude that the equal-time 
Wigner formalism in the limit of a time-dependent electric field E(i) = E(t)e3 is 
totally equivalent to the corresponding quantum Vlasov equation. Most notably, 
the one-particle distribution function can be expressed in terms of the equal-time 
Wigner components:^^ 

m [fB(p, t) - fBvac(p)] + P ■ [v(P, t) - Vvac(p)] 

^(P, t) = . (4.67) 

Alternatively, it is possible to express the equal-time Wigner components w(p, t) in 
terms of J-{p,t), Q{p,t) and Ti.{p,t) as well. As a matter of fact, it turns out that 
only 7 out of possible 16 equal-time Wigner components w{p,t) take non-vanishing 
values: 



ffl(p,t) 






(4.68) 


vi(p,i) 




ga(p,i)] , 


(4.69) 


V2(p,t) 




ge(p,i)] , 


(4.70) 


V3(p,i) 




-G{p,t)] , 

P3 


(4.71) 


ai(p,t) 


= W(P,i) , 




(4.72) 


a2{p,t) 


= T-(-[P,t) , 

2?7i ^ , , , 
= — 7^(p,t). 




(4.73) 


iti,3(p,i) 




(4.74) 



^^Note again, that any function which is considered on a trajectroy /(q, t) can be expressed 
in phase space /(p, t) by replacing the time-dependent kinetic momentum on a trajectory by the 
phase space kinetic momentum, 7r(q, t) p: 

/(P,*) = /(q:*)lq-P+eA(t) ■ 
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4.3 Analytically solvable electric fields 

It is in fact possible to derive an analytic expression for the one-particle distribution 
function .F(q, t) for some specific spatially homogeneous, time-dependent electric 
fields E(i) = E{t)e^. This in turn means that an analytic expression for the equal- 
time Wigner components w(p, t) can be derived as well due to the fact that the 
equal-time Wigner formalism is equivalent to the quantum Vlasov equation in that 
case. 

It has been shown previously that the one-particle distribution function .F(q, t) 
is defined via the Bogoliubov coefficient /3(q, t): 

:F{q,t) = mq,t)W (4.75) 

with: 

/3(q, t) = -ie^G^^^ (q, t) [| + ^c^(q, t)] (q, t) . (4.76) 

Accordingly, an analytic expression for t) can be derived once the mode func- 
tions g^^^ (q, t) as well as the adiabatic mode functions G*-^-* (q, t) are known. As a 
matter of fact, the fundamental system of the equation of motion of a time-dependent 
harmonic oscillator: 



+ uj\ci,t)+ieE{t) 5(q,t)=0, (4.77) 



consists of the mode functions g^^\q,t). Note, however, that this differential equa- 
tion is only for a restricted number of electric fields exactly solvable, most notably the 
static electric field E{t) = Eq as well as the pulsed electric field E{t) = Eq sech^(:^). 
The adiabatic mode functions G^^\q,t), on the other hand, can be calculated once 
the dynamical phase &{q,to,t) is analytically computable: 

G(±)(q,i) = , ^ ^ ^ with e(q,io,t) = / uj{q,t')dt' . 

V2w(q,t)[a;(q,t) T7r3(g3,t)J Jto 

(4.78) 

Given that the one-particle distribution function t) is derivable according to 
this recipe, the auxiliary function G{<l,t) and 7^(q, t) can be calculated as well: 



so that the equal-time Wigner components w(p, t) result from Eq. (4.68) - (4.74). 
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4.3.1 Static electric field 

The static electric field E{t) = Eq can be represented in terms of tlie vector potential: 



A{t) = -Eot , 

so that the equation of motion for the mode functions reads:^"^ 



£ + 4 + [93 + eEot]^ + ieEo g{q, t) = . 



This equation turns into the parabolic cylinder differential equation 



,24 



£ + lu' + l{l + rj) giu)=0 



when performing the variable transformation 



25 



^ [q3 + eEot] with u 6 (-oo, oo) , 



as well as introducing the following shorthand notations: 



1^ = 1 + 31 

eEo e eEo 



with 



The fundamental system consists of parabolic cylinder functions: 



ue 4 



g(+){u) = iV(+)Z?_ 

'-^ 2 

£,(-)(u) = M-^D_,r,{-ue'-f) 



(4.81) 

(4.82) 

(4.83) 

(4.84) 

(4.85) 

(4.86) 
(4.87) 



with iV(±) being normalization factors. The asymptotic behavior of the mode 
functions for u — > — cxd is given by:^^ 



( + )/ N U^-OO 

T '[u) — > 



Ar(+)ieiK+2'?in(|»|)He'r 

\u\ 



(4.88) 
(4.89) 



^^The dependence on the orthogonal momentum will not be explicitly indicated in the fol- 
lowing as the differential equation depends on it only parametrically. 
^^For details on parabolic cylinder functions see [113], chapter 19. 

^^Note that the mode functions g^^^u) will depend only on u due to the linear relation between 
qs and t. 

^^The normalization factors will be chosen in such a way that the mode functions and the 
adiabatic mode functions coincide for n — > — oo: 

= G'*'(m) for u^-oo. 



''Here, without loss of generality _Eo > so that u — > — oo corresponds to f 
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The adiabatic mode functions, on the other hand, are defined as: 

P=F«0(«o,«) 

G^^\u) = , (4.90) 

' rr?e\Jlr] + u^\\j2r\ + =F v\ 

with the dynamical phase being given by: 

1 



e(wo,^i) = 1^ \ ^v^2t] + u2 . (4.91) 

As mentioned previously, the actual value of does not really matter as it only 
fixes an arbitrary phase at a certain instant of time. Choosing the symmetric point 
uq = 0, for instance, one finds :2^ 



9(0, u) = \ [uy/2r] + u2 + 277 ln(u + -^Tlri+v?) - ?? ln(27?)] . (4.92) 
The asymptotic behavior of the adiabatic mode functions for u — > — oo is given by: 

yj2\u\m'^e 

G^-)[u) 1 ^-|p+2»;ln(l«|)+r,+rjln(|)] ^ ^^^^^ 

Accordingly, the properly normalized mode functions read: 

g^^\u) = ^e^K'''"(iHe-^D_,^^(-.e-x), (4.95) 
5(-)(^) = _i^-i[^+^'<)]e-YD,,i-ue^^. (4.96) 

Consequently, it is possible to derive an analytic expression for T{u) = 2|/3(ti)p in 
terms of parabolic cylinder functions, which is displayed in Fig. 4.1: 



y2j;+u'^^ '- ^ 2 2 



(4.97) 



^^The integral is given by: 

y Au\/2ri + = \u\/2ri + + Tyln(u + \/2rf+v?) +J7ln(2) . 

^^The derivative of the parabolic cylinder functions with respect to u, which appears in the 
definition of (i{u), is given by: 

ll^-i+in ( - ne-"^) = ( - ue""? ) + e'^^ ( - ue-"? ) . 
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_J , , \ , , , , , , , , , , , , , , L Q 

-5 5 10 15 

Figure 4.1: One- particle distribution function J^{u) for e = 1 and q± = 0. Starting off with 
J^{u — oo) = 0, the asymptotic value !F{u —>■ oo) = 2e~^^ is reached after some transient 
oscillations. 



4.3.1.1 Pair creation rate 

It has already been mentioned that the first term in the expression for the vacuum 
decay rate: 

(4.98) 



P vac. = } exp — 

n=l ^ 

can be identified with the pair creation rate J\f[e'^e~]:^'^ 

Af[e+e-] = r Au r At^J^'U) = lim -^^^ H di^J^lu) . (4.99) 

Here, the prime denotes the derivative with respect to u. Accordingly, the asymp- 
totic behavior of the parabolic cylinder functions for u — > oo is needed: 

i^-i+.|(-«e-f) _^e-iK+2''MM)]e-^ , (4.100) 

D,a{-ue-f) eH-^+^'^Kl-Dle-^ , (4.101) 

according to which: 

lim J^{u) = 2e-'"'' . (4.102) 



■^"The pair creation rate, i.e. the change in the number of created particles per volume and time, 
is defined as: 
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Accordingly, one finally obtains for the pair creation rate: 



47^3 



exp 



(4.103) 



In Fig. 4.2, the ratio between the pair creation rate A/"[e^e~] and the vacuum decay 
rate P[vac.] is displayed. The difference between those two is rather minor in the 
sub-critical field strength regime e < 1 whereas it becomes sizeable for super-critical 
field strengths e > 1. 

4.3.1.2 Equal-time Wigner components 

It has been mentioned that knowledge of the one-particle distribution function allows 
for determining analytic expressions for the equal-time Wigner components as well. 
In order to simplify notation, the following shorthand notations are introduced:^"'^ 



Viiu) = \D_ 



1+^ 



VJu) 



\Din(- 

2 



ue 4 



ue 4 
2 



-]D 



ue i ] + c.c. 



e 4 Dir, { — ue i)D 



ue 4 



c.c. 



(4.104) 
(4.105) 
(4.106) 
(4.107) 



so that J-{u) can be written as: 



— ZL2 

e 4 



\Vi{u)+[l + 



\V2{u) 



Mu) 



(4.10J 



^Note that T>4,{u) is not yet required but will be used when calculating Q{u) and 7^(u) 
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Due to the fact that: 



i^V-siu) =-uVi{u) , 



as well as: 



Q{u) 



(4.109) 
(4.110) 
(4.111) 

(4.112) 



277 + m2 ' 

it is possible to calculate Q(u) and Iriiu) according to Eq. (4.79) and Eq. (4.80), 
respectively: 



-Hiu) 



2ri -1^L 



(4.113) 



(2r,+«^) 



-W( 1 - e 4 



RV^{u)+V2iu)-^^^^^V,{u) ) . (4.114) 



The transformation back to phase space is required in order to finally calculate the 
equal-time Wigner components as shown in Eq. (4.68) - (4.74). The corresponding 
change of variables q ^ p + eA{t) reads: 



and qj_ ^ P_L , 



so that:^^ 



(4.115) 

I^3(P)], (4.116) 
(4.117) 
(4.118) 



Accordingly, all equal-time Wigner components only depend on the phase space 
kinetic momentum p but not on the time variable t. This is due to the fact that 
the electric field E{t) = Eq is static as well. 



■^^Note that I'i(p) are defined in such a way that they are obtained from 'Di{u) according to: 
I>,(p) = Vi (u J ^Ps) with q± ^p± ■ 
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4.3.2 Pulsed electric field 

The pulsed electric field E{t) = Eosech^^^) can be represented in terms of the 
vector potential: 

A{t) = -£;oTtanh(i) , (4.119) 
so that the equation of motion for the mode functions reads: 

^ + 4 + + e£:oTtanh(^)]2 + ieEosech^^)] g{q,t) = . (4.120) 

In order to actually solve this equation it seems advantageous to introduce the 
dimensionless time variable :^^ 

w= i[l + tanh(^)] with u e [0,1] , (4.121) 

so that: 34 

■^3{q3,u) = q3 + eEoT{2u-l) , (4.122) 

uj\q,u) = el+7Tl{qs,u) . (4.123) 

Accordingly, the equation of motion can be written as: 

[^n(l - u)£[u{l -u)£]+ Lu\(i, u) + UeEQu{l - u)] ^(q, w) = . (4.124) 

Introducing an ansatz for the mode function: 

ff(q,w) =M-'^^(l--«)'^^/i(q,-u) , (4.125) 

the equation of motion turns into the hypergeometric differential equation:^^ 

[l-u)£,+ (c(q) - [a(q) + 6(q) + l\u)^- a(q)6(q)] /i(q, u) = , (4.126) 

with u-independent parameters: 

a(q) = - ieEQT^ - '^^^ + '^^^ , (4.127) 

6(q) = 1 + ieEoT^ - + , (4.128) 

c(q) = l-iTa;(q,0) . (4.129) 



^^The asymptotic limit u — > 0^ corresponds to t — > — oo whereas u — > 1 corresponds to t ^ oo. 
^^Note that the dependence on the orthogonal momentum qx will not be indicated again. 
'^'''For details on hypergeometric functions see [113], chapter 15. 



4.3. Analytically solvable electric fields 



53 



The hypergeometric differential equation is solved in terms of hypergeometric func- 
tions so that the mode functions read:"^^ 

g'-+\q,u) = N^+\q)u-'-^^^{l-u)^^^^F{a,b,c;u) , (4.130) 
g^-\q,u) = 7V(-)(q)u^^^^(l - w)-^^^^F(l - a,l - 6,2 - c;u) . (4.131) 

Accordingly, the asymptotic behavior of the mode functions for u ^ 0^ is given by: 

g(+)(q,„) Ar(+)(q)e-^M«) , (4.132) 

g(-){ci,u) 7V(-)(q)e^i°(") . (4.133) 

In analogy to the static electric field E{t) = Eq, the asymptotic behavior of the 
adiabatic mode functions: 

q=ie(q,Mo,«) 

G(±)(q,n) = , , (4.134) 

^/2uJ{q,u)[uJ{q,u) iF7r3(q3,w)] 

with: 

e(q,uo,u) = l- rdz;^^^^, (4.135) 



is determined in order to fix the value of the normalization constants A^'-^^(q). The 
dynamical phase can in fact be split into a logarithmically divergent part as well as 
an irrelevant regular part $(q, mo,0) in the asymptotic limit u — 0^:'^^ 

e(q, uo,u) ^^^^ Hu) + $(q, uo, 0) , (4.136) 

so that the asymptotic behavior of the adiabatic mode functions for u ^ 0^ is given 
by: 

G(±)(q,z,) ^ e^^-M ,^^Hu) (4^137) 



^^The fundamental system of the hypergeometric differential equation {/i'^' (q, u), h^- \q, u)} in 
the neighborhood of u = 0'*' consists of: 

(q, u) = F{a, b, c; u) and /i*"' (q, u) = (1 _ n)-'™<'^'''F(l - a,l - b,2 - c;u) . 

■^^The integral is given by: 
J du^^^ = -i7r[(.j(q,0)+a;(q,l)] - 2e£;or [ln(2) + In (^7r3{q, n) + Lj(q, u))] 

+ a;(q, 0) [ln(u) + In (^^^^^ j - In (^ei + 7r3(q, 0)7r3(q, u) + a;(q, 0)a;(q, u)^j 

- a;(q, 1) [ln(l - n) + In (^^i^) - In + 7r3(q, l)7r3(q,J.) + a;(q, l)a;(q, u))] 
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Figure 4.3: Time evolution of the one-particle distribution function !F{q,u) for e = 1, 
T = — and qx = 0. Starting off with .F(q, u —>■ 0+) = 0, the asymptotic value .F(q, m ^ 1) 
is finally obtained. Note that all the remaining figures in this section are displayed in terms 
of the phase space kinetic momentum p-^ instead of the canonical momentum q^. 

The normalization constants A'^'^^-'(q) are accordingly given by: 

P=Fi*(q,M0,0) 

iV(±)(q)= (4.138) 
v/2^(q,0)[c^(q, 0)^^3(93,0)] 

Consequently, it is possible to derive an analytic expression for J-{ci, u) = 2|/3(q, w)p 
in terms of hypergeometric functions, which is displayed in Fig. 4.3: 

^(q,^) = |Ar(+)(q)|2[i + !gg^]|^i(q,^)+i^2(q,«)|', (4.139) 

with: 

.Fi(q,u) = 1^(1-^)^^(1 + 0,1 + 6,1 + 0;^) , (4.140) 

J-2(q,w) = [a;(q,u)-(l-u)w(q,0)-wu;(q,l)]F(a,6,c;n) . (4.141) 



4.3.2.1 Asymptotic one-particle distribution function 

The analytic expression for one-particle distribution function in terms of hyperge- 
ometric functions drastically simplifies in the asymptotic limit. The asymptotic 

^^The derivative of the hypergeometric function with respect to the last argument is given by: 

■§^F{a, b, c; u) = ^F{1 + a, I + h,l + c;u) . 
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Figure 4.4: Asymptotic one-particle distribution function JF(q, u ^ 1) for e = 1, r = — 
and qx = 0. Note again that the interpretation of J^(q, u) as momentum distribution of 
real particles is in fact only possible at asymptotic times u I. 



behavior of the hypergeometric functions for u — > 1 is in fact given by:^^ 

F{a,b,c;u) "-=1 , (4.142) 

Fil+a,l+b,l+c;u) (i - , (4.143) 

so that: 

2 



Im J'^(q, u) - . (4.144) 



2b r(i+c)r(i+a+b-c) 



c r(i+a)r(i+fe) 



Applying the transformation formula for gamma functions, one finally obtains a 
simple analytic expression for the asymptotic one-particle distribution function:^'' 

2sinh(^[2e£;oT+(.^(q,0)-a;(q,l)])sinh(^[2e_Eor-a;(q,0)+a;(q,l)]) r a i 
Sinn I 7rra;(q,U) ) smh I 7rTa;(q,l) I 

In Fig. 4.4, the asymptotic one-particle distribution function J-{q,u — > 1) is dis- 
played which is obtained once the electric field vanishes. 



^^Note tiiat tiie liypergeometric functions have been transformed by means of the linear trans- 
formation formula: 

F(a,f3,r,u) = Eh)Ih^^F(a,f],a + f]-^ + 1-1 -u) + 

(1 _ „)7--/3rMnfd|p)F(^ -a.rf-Prf-a-P + l;l-u) . 

''"For details on gamma functions see [113], chapter 6. The following relations are used here: 
r(l + a)=ar(a) and |r(l + = ^^^^ . 
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T(q,u^l)vs. 21n[l-'W(q)] 




4.3.2.2 Particle density vs. imaginary part of the eflPective action 

It is very instructive to compare the resuh for the density of created particles 
A/'[e^e~], which is obtained from the quantum Vlasov equation, with the result 
for the imaginary part of the effective action 2ImiSefr, which has been obtained 
from semi-classical methods previously [18].^^ The density of created particles is 
given by: 



M[e+e-] = lim / — i-^(q,^x) , (4.146) 
whereas the imaginary part of the effective action density reads: 



with: 



2Im5eff = J ^21n[l - W(q)] , (4.147) 

^ sinh (^[2ei;or+t^(q,0)-a;(q,l)]) sinh ( ^ [2ei;or-a;(q,0)+^(q.l)]) 
sinh (^[2e£;or+w(q,0)+w(q,l)]) sinh (^[2e_Bor-w(q,0)-w(q,l)]) 

In Fig. 4.5, the integrand of Eq. (4.146) is compared with the integrand ofEq. (4.147). 
One finds that 21n[l — W(q)], which is related to vacuum decay probability, is in 
fact always larger than the asymptotic one-particle distribution function J-{u — 1), 



In contrast to the static electric field E{t) = Eq where it has been only sensible to talk about 
the pair creation rate J\f[e^e~], i.e. the change in the number of created particles per volume and 
time, it is now appropriate to consider the density of created particles A/'[e^e~], i.e. the number of 
particles per volume, which are created during an electric pulse of finite duration: 



M[e+e-] = dt I j^T(ci,t) = lim I 



-JP(q,t) 
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which characterizes the momentum distribution of created particles. This is in ac- 
cordance with the discussion for the static electric field E(t) = Eq where it has 
been pointed out that the vacuum decay rate P[vac.] is always larger than the pair 
creation rate A/'[e^e~]. 

In Fig. 4.6, the ratio between the density of created particles A/"[e~'~e~] and the 
imaginary part of the effective action density 2Im<Seff is displayed. Similarly to 
the static electric field, cf. Fig. 4.2, the difference between those two quantities is 
again rather minor in the sub-critical field strength regime e < 1 whereas a sizeable 
deviation occurs for super-critical fields e > 1. 

4.3.2.3 Equal-time Wigner components 

Given the analytic expression for the one-particle distribution function in Eq. (4.139), 
it is again possible to derive analytic expressions for the equal-time Wigner compo- 
nents as well. As a matter of fact, one has to calculate Q{q,u) and Ti.{q,u) first, 
which are given by: 



^(q^") = ^y^^'(q.^)> (4-149) 

W(q, u) 



2eEoe_\_u(l — u) 



1 - ^(q,n) - -'^'^lllf'-^ :F'{ci^u) + [ig^] V-(q,.)] . (4.150) 

Here, the prime denotes differentiation with respect to u again. Note, however, that 
both !F'{q^,u) and !F"{c\,u) yield rather cumbersome expressions as there occurs no 
simplification like for the static electric field E{t) = Eq. 
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The change of variables q — > p + eA{t) then allows for the transformation back to 
phase space: 

q3^P3-eEoT{2u-l) and q± ^ P_l , (4.151) 

so that: 

■Ks{q,u)^P3 and w(q, u) ^ u;(p) . (4.152) 

Accordingly, any function which depends only on the canonical momentum qs shows 
a dependence on both the phase space kinetic momentum and the time variable 
u after this transformation:^^ 

t^(q,0) ^ ^ei + [p3-2ue^or]2 , (4.153) 
u;(q,l) ^ ^el + [p3 + 2{l-u)eEoTP . (4.154) 

4.4 Pulsed electric field with sub-cycle structure 

In the previous section, an analytic expression for the one-particle distribution func- 
tion J-{ci, t) in the presence of a static electric field E{t) = Eq as well as a pulsed 
electric field E{t) = i?osech^(|:) has been derived. For an arbitrary electric field 
E{t), however, no such analytic solution can be found so that the investigation of 
the Schwinger effect has to be based on numerical results.''^ 

As a matter of fact, both the static and the pulsed electric field are not very 
realistic in the sense of representing a field configuration to be realized at high- 
intensity laser facilities: Due to the fact that the field strengths which are required 
to create electron-positrons pairs out of the vacuum are most probably produced in 
the focus of colliding laser pulses, the electric field shows in fact both spatial and 
temporal variations. '^^ However, as the spatial focussing scale is much larger than 
the Compton wavelength Ac, the effect of spatial variations might be ignored in first 
approximation, resulting in an electric field configuration which is well suited for an 
investigation by means of the quantum Vlasov equation. 

^^Most notably, the parameters a(q), b(q) and c(q) appearing in the first three arguments of the 
hypergeometric functions depend on both pa and u after this transformation. The transformation 
of the dimensionless time variable u to the usual time variable t is obtained via Eq. (4.121). 

''■^A collection of Fortran solvers for the initial value problem for ordinary differential equation 
systems named ODEPACK, which has been invented by the Lawrence Livermore National Labo- 
ratory, has been applied in this work, cf.: http://www.netlib.org/odepack/ 

^''Magnetic fields are neglected as I assume two counter-propagating laser pulses, forming a 
standing-wave electric field in such a way that the magnetic field vanishes. 
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Figure 4.7: Shape of the electric field for = and different values of (t = 3 (red), a = A 
(purple) and a = 5 (blue). The Gaussian envelope is given as reference (dotted black). 

A simple model of such an electric field in the focus of two colliding laser pulses, 
which takes into account both the finite pulse duration as well as the sub-cycle 
structure of a laser pulse, is given by: 



Here, u denotes the laser frequency, t defines the total pulse length and (p is the car- 
rier phase which accounts for a possible shift between the maximum of the Gaussian 
envelope and the maximum of the sub-cycle oscillation.'^^ Note that the appearance 
of a total of three scales makes the investigation of the Schwinger effect in such a 
type of electric field rather involved. It is convenient to introduce the parameter: 



which is a measure of the number of oscillation cycles within the Gaussian envelope. 
The corresponding vector potential can be expressed in terms of error functions:^^ 



The electric field Eq. (4.155) is displayed for 99 = and different values of a in 
Fig. 4.7. Additionally, it is shown for a = 5 and different values of ip in Fig. 4.8. 

^^These studies have in fact been strongly motivated by the sensitive carrier phase dependence 
of strong- field ionization experiments in atomic, molecular and optical physics [114]. 

*^The electric field strength Eo, the laser frequency uj as well as the pulse length t have to be 
considered with respect to the intrinsic scale of QED, which is set by the electron mass m. 

^^For details on the error function see [113], chapter 7. 




(4.155) 



a = UJT , 



(4.156) 




(4.157) 
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4.4.1 Momentum distribution 

As a first issue the momentum distribution of created particles is investigated. To 
this end, the quantum Vlasov equation in its differential form Eq. (4.36) - (4.38) 
is numerically solved with the electric field and vector potential being given by 
Eq. (4.155) and Eq. (4.157), respectively. Note that the backreaction mechanism is 
neglected as the focus lies on the subcritical field strength regime e < 1. 

The parameters are actually chosen in a such a way that t = corresponds 
to a total pulse length of several times 10~^^ s which lies in the anticipated range of 
future XFELs or may become realizable with higher harmonics or secondary beam 
generation of optical lasers. Moreover, the field strength parameter e = 0.1 and the 
orthogonal momentum qj^ = 0. 



4.4.1.1 Carrier phase (p = 

The Schwinger effect in the presence of the electric field E{t) with vanishing carrier 
phase If = has been investigated by means of a Gaussian WKB approximation 
previously [34]. It has been shown that the asymptotic momentum distribution in 
the non-perturbative regime is given by:"^^ 

^WKB(q)~exp(-f[l-f]-^[f(?| + qi]) , (4.158) 

^^The non-perturbative regime is characterized by a Keldysh parameter 7-^1. Note, however, 
that the appearance of two different time scales i and t makes the definition of 7 ambiguous [35]. 
Considering the laser frequency u! as dominant time scale, which is justified at least for a > I, 
suggests the definition: 
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Figure 4.9: Asymptotic one-particle distribution function J^{q,t cxd) when passing from 
CT = 3 (red) over a = A (purple) to cr = 5 (blue). The main peak shifts to smaller kinetic 
momenta and distinct oscillations in tune with the laser frequency uj set in when the 
number of oscillation cycles is increased. All other parameters are given in the text. 

with 7^ = (1 + ^)7^- As it turns out, .FwKB(q) misses in fact several charac- 
teristic properties of the momentum distribution. Before actually discussing these 
shortcomings of the Gaussian WKB approximation, I focus on the asymptotic one- 
particle distribution function JF(q, t —>■ oo) which is displayed in Fig. 4.9. It is 
obvious that various effect occur upon passing from cr = 3 to cr = 5 which corre- 
sponds to increasing the number of oscillation cycles within the Gaussian envelope: 

First of all, the main peak of the momentum distribution shifts to smaller 
kinetic momenta and becomes narrower upon increasing the value of cr: As a mat- 
ter of fact, the electric field Eq. (4.155) behaves much like a single pulse for small 
values of a so that all created particles are accelerated into one direction only. By 
increasing the value of cr, however, the net acceleration of created particles by the 
electric field becomes zero, resulting in a momentum distribution which is peaked 
around: 



The momentum distribution shows in addition distinctive oscillations, with the 
oscillation scale set by the laser frequency uj. Most notably, these oscillations become 
even more pronounced for increasing values of cr. It has to be emphasized that this 
striking feature is totally missed by the Gaussian WKB approximation as displayed 
in Fig. 4.10. Additionally note that the width of the momentum distribution is 
predicted somewhat broader by the Gaussian WKB approximation compared to 




(4.159) 



.F(q, t oo). 
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In order to qualitatively understand the oscillatory behavior of the momentum 
distribution, it is convenient to recall that the Schwinger effect in a spatially homo- 
geneous, time-dependent electric field E{t) can be reformulated as a one-dimensional 
quantum mechanical over-barrier scattering problem, with the scattering potential 
given by — a;^(q, [17, 19, 115]. According to this approach, the reflection coef- 
ficient |7?-(q)p of the scattering problem is related to the asymptotic one-particle 
distribution function J-(q^,t — > oo): 

Accordingly, the asymptotic one-particle distribution function changes J-{q, t — s- oo) 
as the shape of the potential — a;^(q, t) varies for different values of q. Therefore, the 
oscillatory behavior of the momentum distribution can be interpreted as resonance 
phenomenon in the equivalent one-dimensional scattering problem. As a matter of 
fact, this also explains why the spacing between the local maxima is in tune with 
the laser frequency oj. 



As the Gaussian WKB approximation misses some characteristic features of the 
asymptotic one-particle distribution function J-{ci,t oo), attempts have been 
made to go beyond this approximation. One possible way is to use the WKB in- 
stanton action approach, according to which the momentum distribution of created 
particles is given by [21, 41, 43]: 

•^WKB(q) - . (4.161) 
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Here, iSinst.(Q) denotes the instanton action. It has been shown that the leading 
order instanton action can be defined in the complex t-plane as a contour integral:^^ 

2Sinst.{ci)=i(^^dtuj{(i,t)=ij^^dt^el + [q3-eA{t)]'^ , (4.162) 

with the path T being chosen around the branch cut in the complex t-plane. The 
branch points t± are defined according to:^*' 



(4.163) 



Performing a change of variables 



e 2 rErf 



+ 



V2t ^2 



+ c.c. 



(4.164) 



the instanton action reads: 



25inst.(q) 



dT ; — exp — K 

cos{ujt) ^\2t^ 



(4.165) 



with t being considered as function of T.^^ This expression exhibits a number of 
singularities: The exponential function as well as the square root diverge for large 
T whereas the cosine term exhibits poles which are determined by; 



(2n - l)7r 
2u! 



(4.166) 



Expanding the integrand of Eq. (4.165) in a Laurent series: 



52 



oo 

-l+l \ ^ 



1=0 k=0 



u;'.,„,2 5a; 4 610;^ 



720 



t{Tf 



, (4.167) 



''^Actually, Eq. (4.162) gives the instanton action for sQED. It has been shown, however, that 
the instanton action for sQED can also be apphed in QED [41]. 

^''For a static electric field with A{t) = —Eot this equation has two unique solutions. For the 
vector potential Eq. (4.157), however, this equation becomes ambiguous. Note that this ambiguity 
in the definition of the branch points will be disregarded in the following. 

In fact, Eq. (4.164) cannot be inverted analytically but only by reversion of the corresponding 
series. The first few terms in this expansion are given by: 

*(T)^T+i±^T3+I±iig^T^ + 0(TV 



6t2 



120t4 



The detailed calculation can be found in Appendix A. 6. 

Note that only the poles for n = and n = 1 are taken into account in this expansion whereas 
the contributions from all the other n are disregarded. As it turns out, the omission of those 
additional contributions results in a momentum distribution which misses the oscillatory behavior 
at the end. 
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the instanton action can be calculated by summing up the existing residues. Here, 
the coefficients C;(q) are determined in terms of the series expansion: 



^^-^p -2^G(q)T . (4.16J 



1=0 

Introducing the following notation for the instanton action: 

25,n..(q) = ^Ef^' (4-169) 
the first terms are found to be given by: 

5W(q) = 1 , (4.170) 
5{^>(q) = i^(5i-K) > (4-171) 
5W(q) = fe^ - i^iei + M) , (4.172) 

In Fig. 4.11, the asymptotic one-particle distribution function .F(q, 00) is compared 
with the Gaussian WKB approximation .FwKB(q) and the improved WKB approx- 
imation .^■v\rKB(l) cr = 5. The improved WKB approximation fits very well the 
averaged momentum distribution. However, neither the Gaussian nor the improved 
WKB approximation predict the oscillatory behavior. It was only very recently that 
this discrepancy has been explained in detail by taking into account interference ef- 
fects in an extended WKB approximation [38, 40]. 



4.4. Pulsed electric field with sub-cycle structure 



65 




Figure 4.12: Asymptotic one-particle distribution function .F(q, t —>■ oo) for tp = j (top) 
and </' = f (bottom). All other parameters are given in the text. 

4.4.1.2 Carrier phase f ^ 

I now draw my attention to the Schwinger effect in the presence of the electric field 
Eq. (4.155) with non-vanishing carrier phase ip ^ 0. In fact, this field configuration 
has not been considered before as a non- vanishing carrier phase breaks the symmetry 
E{t) = E{—t) of the electric field, which in turn makes the imaginary time treatment 
of the WKB approximation significantly more complicated. 

In Fig. 4.12, the asymptotic one-particle distribution function J-{ci,t oo) is 
displayed for different values of ip.^^ Apparently, an increasing phase offset makes 
the oscillatory behavior of the momentum distribution more pronounced. The most 
distinctive signatures are in fact found for tf = ^ when JF(q, t — > oo) actually van- 
ishes at the minima of the oscillations. Moreover, the main peak shifts to momenta 

7^ again, even for large values of a. 

^^Equivalently, a negative value of ip < results in a momentum distribution which is mirrored 
at p3 = 0. 
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This characteristic behavior of the asymptotic one-particle distribution function 
.F(q, t — > Qo) can be explained in a simple way by considering the equivalent one- 
dimensional scattering problem again: The electric field is symmetric with respect 
to time reversal for ip = whereas it becomes antisymmetric for (p = ^. Note that 
the corresponding vector potential A{t) shows exactly the opposite behavior: 

y,=0: E(t) = E{-t) and A{t) =-A{-t) , (4.174) 
(/3=f: E{t) =-E{-t) and A{t) = A{-t) . (4.175) 

Owing to the symmetry property of A{t), the scattering potential is in fact sym- 
metric with respect to time reversal for = ^: 

(/.= §: u;2(q,f) =a;2(q,-t) . (4.176) 

As such symmetric scattering potentials allow for perfectly transmitted states :^^ 

|7^(q)|2 = , (4.177) 

it is possible to trace the distinctive signatures in the momentum distribution back 
to a resonance phenomenon: The momentum values at which t oo) vanishes 
are the same at which perfectly transmitted states in the equivalent one-dimensional 
scattering problem occur. Note that a more detailed explanation of the oscillatory 
behavior has recently been given by means of an extended WKB approximation 
[38, 40]. 

In analogy to strong-field ionization experiments [114, 116, 119, 120], these dis- 
tinctive signatures in the momentum distribution may serve as sensitive probe of 
sub-cycle structure in upcoming high-intensity laser experiments. In addition to the 
density of created particles, these results suggest a number of new observables such 
as the main peak position jj^^ , the width of the momentum distribution and, most 
notably, its oscillatory structure. In particular, the characteristics of these oscil- 
lations provide precise information about the carrier phase ip and the total pulse 
length parameter t. As the latter features are hard to control a priori in an abso- 
lute manner for a high-intensity laser, these momentum signatures could serve as a 
tomograph of the laser pulse, providing for a unique means to verify and confirm 
design goals of future high-intensity laser systems. 

^*This argumentation is guided by intuition from simple one-dimensional scattering problems 
such as the square well. Unluckily, I am not aware of any rigorous proof of this statement. 

^^Note that there is a close similarity between this effect and the matterless double slit experiment 
[116, 117, 118]. 
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Figure 4.13: Comparison of .Fsp(q, i —>■ oo) (blue) with .Fsc(q, t —>■ oo) (purple) for ip = 
{top) and V = f (bottom). The momentum distributions are only shown for kinetic 
momenta > pg*^. All other parameters are given in the text. 

4.4.2 Quantum statistics effect 

So far. all numerical simulations have been based upon the quantum Vlasov equation 
in the framework of QED:^^ 

t 

•^sp(q,i) = 5Qsp(q,i) J Qsp(q, " 2^sp(q, 0] cos [2e(q, i', t)] . (4.178) 

^vac 

*^Note, that the quantum Vlasov equation has now been written in terms of the one-particle 
distribution function for one spin direction jFsp(q, t) in order to recognize the close similarity 
between QED and sQED later on. Note that Eq. (4.33) is recovered according to: 
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The analogous quantum Vlasov equation in the framework of sQED reads [29]: 

t 

■^sc(q,i) = 5Qsc(q,0 j Qsc(q,i')[l + 2^.c(q,0]cos [2e(q,i',0] , (4.179) 

^VUC 

with: 

Qsp(q,0 = ^ and Q.(q,0 = ^^^t^. (4.180) 

There are two major differences between the quantum Vlasov equations in the frame- 
work of QED and sQED, respectively: First of all, a Bose enhancement factor ap- 
pears for scalars [1 + 2.7-'sc] whereas a Pauli-blocking factor is found for spinors 
[1 — 2J-sp\- Moreover, the different prefactors Qs-p{<\, t) and (5sc(q5 1) account for the 
difference in the phase space occupation between fermions and bosons. 

In Fig. 4.13, the asymptotic one-particle distribution function .7-'sp(q, t — > oo) for 
spinor particles is compared with the asymptotic one-particle distribution function 
.7-'sc(q, t — > oo) for scalar particles. The averaged momentum distribution is in fact 
the same for both QED and sQED, however, quantum statistics plays a crucial rule 
for the detailed shape: The momentum distribution for spinor particles shows a 
local maximum at those momentum values at which the momentum distribution for 
scalar particles shows a local minimum, and vice versa. This out-of-pliase behavior 
between QED and sQED is in fact found for any value of ip. As for the explanation of 
the oscillatory structure of the momentum distribution, this out-of-phase behavior 
between QED and sQED has very recently been explained in detail in an extended 
WKB approximation [38, 40[. 



Chapter 5 

Schwinger effect for E(x,t) 



In this chapter I present the results of an ab initio simulation of the Schwinger 
effect in the presence of a simple space- and time-dependent electric field. This 
investigation represents in fact substantial progress as there have not been any 
rigorous studies of the Schwinger eff"ect in space- and time-dependent electric fields 
so far. Note, however, that a solution of the equations of motion for the equal-time 
Wigner components Eq. (3.36) - (3.43) in the framework of 3 + 1 dimensional QED 
(QED3^]^) is currently not feasible for computational reasons. Accordingly, I have 
to restrict myself on 1 -|- 1 dimensional QED (QED^^^) within this thesis.^ 

In Section 5.1, I briefly introduce the equal-time Wigner formalism in QED]^_|_]^ 
upon which the subsequent simulations are based. In Section 5.2, various strate- 
gies for solving the equations of motion for the equal-time Wigner components are 
presented, including both an exact and an approximate scheme as well as a local 
density approximation. Eventually, in Section 5.3, I discuss the numerical results 
concerning the momentum distribution, the charge distribution as well as the total 
number of created particles. Moreover, a comparison between the exact solution 
and the approximate ones is drawn. 



5.1 Equal-time Wigner formalism in QED^,^^ 

It has been stated that the equal-time Wigner formalism in mean field approximation 
is complete in the sense that knowledge of the equal-time Wigner components suffices 
to calculate any physical observable. Accordingly, the investigation of the Schwinger 
effect can be considered as completed from a theoretical point of view. 

^In order to roughly estimate the required memory for a direct numerical simulation in QEDg.,.^, 
consider a discretization of the PDE system Eq. (3.36) - (3.43) on a finite grid, with the number 
of grid points in each direction being of the order of rii ~ 0(10^). Accordingly, in order to fully 
characterize the state of the system at a given instant of time, the required memory is of the order 
of: 

m ~ 0(10^^6) ~ C>(10TB) . 

Actually, this estimate is far too conservative: The experience in simulating vacuum pair creation 
in 1 + 1 dimensional systems shows that the number of grid points has to be chosen to be m = 
0(10^) - 0(10'*) in order to yield both high accuracy and good convergence. Accordingly, the 
required memory is far beyond the scope of current technology. 



5.1. Equal-time Wigner formalism in QED 



70 



Nonetheless, the numerical simulation of the PDE system Eq. (3.36) - Eq. (3.43) 
poses still a serious problem from a practical point of view: This PDE system is for- 
mulated in terms of 16 equal-time Wigner components "w(x, p, t) which themselves 
depend on six phase space variables {x, p} as well as the time variable t. Accord- 
ingly, any numerical simulation is doomed to fail for computational reasons for the 
time being. 

As a matter of fact, one could restrict oneself to a highly symmetric electro- 
magnetic field F^'^{x) in order to reduce the number of phase space variables by a 
proper choice of the coordinate system. Choosing a cylindrical symmetric configura- 
tion, for instance, the polar angle could most probably be eliminated for symmetry 
reasons.^ Nonetheless, four phase space variables corresponding to the directions 
parallel and perpendicular {x_\_,p_\_) to the symmetry axis, respectively, 

would still remain. 

In order to avoid all these difficulties I restrict myself to QEDj^^j^ from the very 
beginning.'^ Accordingly, one is left with only two phase space variables {x,p\. 
Moreover, the number of equal-time Wigner components ^{x,p,t) can be reduced 
from 16 to 4 as well.^ Note, however, that the notion of a magnetic field does not 
exist in QED]^^]^ in contrast to QED3^]^ so that fermions have to be considered as 
spinless particles. Nonetheless, this peculiarity of QED;^^^ should not be of severe 
consequence: As the focus of this thesis lies on vacuum pair creation in electric 
fields, magnetic effects would have been neglected anyway. Accordingly, even if the 
results from QED]^_|_]^ cannot be transferred one-to-one to QED3_|_]^, it can still be 
hoped that they contain at least the main features of the Schwinger effect in the 
presence of space- and time-dependent electric fields in QED3^]^. 

Moreover, it has to be emphasized that the experience which is gained by the nu- 
merical simulation of the Schwinger effect in QED]^_|_]^ can most probably be carried 
over to the investigation of the Schwinger effect in QED3_|_]^ as well. 

^The generalization of the equal-time Wigner formalism to curvilinear coordinates has not been 
worked out in this thesis. This extension of the current work is surely worthwhile to consider in 
future investigations. 

^For massless fermions this is called the Schwinger model, which can be solved analytically 
[121, 122]. For details on the massive Schwinger model see [123, 124]. 

''Accordingly, the required memory for a direct numerical simulation on a grid with n, ~ 0(10'') 
grid points in each direction is feasible with current technology: 

m ~ C»(10®B) ~ O(IGB) . 
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5.1.1 QED in 1 + 1 dimensions 

The Lagrangian of QED]^_|_]^ looks alike the Lagrangian of QED3_|_]^ given in Eq. (3.1): 

^, A) = i (^^^^[idf, - eA^]* - [idf, + eA^]^j''^^ - m** - {F'^'F^, . (5.1) 

Note, however, that space-time is only two-dimensional so that fi = {0, 1}. Accord- 
ingly, the Dirac algebra is composed of two Dirac gamma matrices only:^ 

{7^ Y} = ^g^"" with [7°]^ = 7° , [7^ = "7^ • (5.2) 

As a consequence, the Dirac field ^'(x) and the adjoint field ^'(x) posses only two 
components. Moreover, it is again possible to define a chirality matrix:^ 

7^ = 7O7I with [7^]^ = 75 . (5.3) 

Additionally, the field strength tensor F'^'^{x) possesses only one non-trivial compo- 
nent which can be identified with the electric field:^ 

F^^(x) = 5^A^(x) - d^Af^ix) with Fi°(x) = ^(x) . (5.4) 

5.1.2 Covariant Wigner formalism 

The derivation of the equation of motion for the equal-time Wigner function in 
QED]^_I_]^ is along the same lines as the corresponding derivation in the framework 
of QED3_|_]^, which has been discussed in detail in Chapter 3. Starting off with the 
covariant Wigner operator: 

yV^^\x, ?) = \j d2ye^P-ye*^^-V^^«^('<+«^)-^ [^(x - ^), *(x + |)] , (5.5) 

the covariant Wigner function is again defined as the vacuum expectation value: 

W'^'^\x,p) = {n\W^'^\x,^)\Q) . (5.6) 

^The Dirac algebra in two dimensions can be represented in terms of the first two Pauli matrices: 

7" = (Ti and = —ia2 ■ 

^With the previous choice of Dirac gamma matrices, the chirality matrix is given by the third 
Pauli matrix: 

7"'^ = 0"3 . 



^In this chapter, sans serif variables denote two- vectors, for example p = = (p^jp) 
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The equation of motion for the covariant Wigner function is then again found by 
taking advantage of the Dirac equation and its adjoint. As the BBGKY hierarchy 
is again truncated at the one-body level in mean field approximation, one finally 
obtains the self-consistent PDE system: 

iD^[7^' _in^{^M^yv;(2)} = -2imW^^^ , (5.7) 

iD47^>V(')}-^^^ [7^ = , (5.8) 

and 

d-^F'^'^ix) = (Oil [*(x),7'^*(x)] |Q) = tr[YW^'\x,p)] , 

with: 

1 

D^{x, p) ^ d^^-e f\d^ F^,{x - i^dp)d; , 

1 

n^(x, p) = - ie j\ diiF^,{x - iid^)d'; . 

•'~2 

Decomposing the covariant Wigner function into its Dirac bilinears: 

W(2)(x,p) = i [S + i7^P + 7''Vj , (5.12) 

one may derive equations of motion for the covariant Wigner components W(x, p). 
Adopting both the commutator and the anticommutator between 7^ and the Dirac 
bilinears:^ 





{7^■} 


[7^■] 


I 


27^* 





75 





-2e'^«7o 




2g^- 


2e'''^75 



(5.13) 



the resulting equations of motion split again into a set of inhomogeneous equations: 

n^V^ = mS , (5.14) 
n'^S - le^'^D^V = mV^' , (5.15) 
D^V - D-'V^ = -2me>"'¥ , (5.16) 

^Again, the usual convention e'^^ = 1 for the totally antisymmetric tensor is used. 



(5.9) 

(5.10) 
(5.11) 
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as well as into a set of homogenous equations: 

Df'Y^ = 0, (5.17) 

D^'S + 2e'"'U^¥ = 0, (5.18) 

n^v^ - n'^v^ = . (5.19) 

5.1.3 Equal-time Wigner formalism 

The equal-time Wigner function is then again defined as the energy average of the 
covariant Wigner function: 

wW(x,p,t)^ J ^W(2)(x,p) . (5.20) 

Decomposing the equal-time Wigner function into its Dirac bilinears: 

(x, t) = i [s + *75P + 7^v,,] , (5.21) 

the equations of motion for the equal-time Wigner components w{x,p,t) are again 
found by taking the energy average of Eq. (5.14) - (5.19). The time- evolution 
equations are then given by:^ 

DifB - 2pp = , (5.22) 

Avo + £v = , (5.23) 

Dtv + £vo = -2m p , (5.24) 

Dtp + 2ps = 2m V , (5.25) 



with 



1 

Dt{x,p,t)= §^+e j\_diE{x + ii§^,t)§^ . (5.26) 



2 

The constraint equations, on the other hand, read 



ffi^^J + ntffi - i^p = mi^o , (5.27) 

Vg^' +IItVo —p^ = ms , (5.28) 

+ V -pvo = , (5.29) 

pW + Htp + if ffi = , (5.30) 



^The pseudoscalar equal-time Wigner component p(x,p, t) in QEDj_,_j takes on the role of the 
tensorial equal-time Wigner component ti(x, p,£) in QED3_,_j as cr^" = 17^ in 1 -|- 1 dimensions. 
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with 

1 

Utix,p,t)=ie l\diiE{x + i^^^,t)-l-^ . (5.31) 

It can be shown again that the transport equations for all higher energy moments 
w["l(x,p, t) are identically fulfilled in mean field approximation once the equal-time 
Wigner components obey Eq. (5.22) - (5.25). Along with appropriate vacuum initial 
conditions: 

Till p 

fBvac(p) = and Vvac(p) = , (5.32) 

u;{p) uj{p) 



with uj{p) = \/m?+p^ , the equal-time Wigner formalism is complete in the sense 
that w(x,p, i) encodes the information of all higher energy moments. 

It has to be emphasized that the equal-time Wigner formalism in QED]^.,.]^ rep- 
resents a drastic simplification compared to the equal-time Wigner formalism in 
QED3^]^ as the number of equal-time Wigner components is reduced from 16 to 4. 
Nonetheless, it can still be hoped that the Schwinger effect in QEDi_^_i contains at 
least the main features of the Schwinger effect in QED3_|_]^. 

5.1.4 Observables and marginal distributions 

The observables in this theory are again found from an analysis of the conservation 
laws. Noether's theorem proves again the conservation of the total charge Q, the 
total energy £ and the total momentum V as well as the Lorentz boost operator 



Q = ej dTwoix, p,t) , (5.33) 

£= J dT[ms,{x,p,t) +p^{x,p,t)] + ^ J dx\E{x)\^ , (5.34) 

V= JdTpMx,P,t), (5.35) 

]C = ~ J dr a; [ms(x,p, i) -h p v(a:,p, t)] — i J dxx\E{x)\'^ , (5.36) 

with dr = dxdp/(27r) being the phase space volume element. These conserved 
quantities suggest to define various marginal distributions, most notably the real 
space charge density q{x,t) and the momentum space charge density q(p,t): 

/dp f dx 

j^vo{x, p,t) and q{p,t) = J j^V(,{x,P,t) , (5.37) 



The notion of angular momentum does not exist in 1 -|- 1 dimensions. 



5.1. Equal-time Wigner formalism in QED 



75 



as well as the real space energy density e{x, t) and the momentum space energy 
density e{p, t) of Dirac particles; 

e{x,t) = I — — - [ms(x,p, + p t)] , (5.38) 



j^[ms{x,p,t) +p^{x,p,t)\ . (5.39) 

The total energy of the Dirac particles Snit) is in fact composed of a vacuum 
contribution as well as a contribution due to created particles. It is hence indicated 
to introduce the total energy of created Dirac particles S^{t): 

£l{t) = Soit) - £d,..c , (5.40) 

with: 

^hi^) = J dT[m[s,{x,p,t) - Svac{p)] +p[^{x,p,t) - ^yac{p)]] ■ (5-41) 

On the other hand, £fj{t) should also be calculable by integrating a particle number 
quasi- distribution n{x,p,t) times the one-particle energy uj{p): 

£Ut) = j dTu;{p)n{x,p,t) . (5.42) 

Correspondingly one defines the real space particle number density n{x,t) and the 
momentum space particle number density n[p,t)'}^ 

f dp m[s(a;,p,i) - Svac(p)] +p[v(a;,p,i) - Vvac(p)] 

" JW) ^) ' ^ ^ ^ 

n(p t) = f m[s(x,p,t) - Svac(p)] +p[v(3:,P,t) - Vvac(p)] 

J (27r) uj{p) 

The total number of created particles is then calculated in terms of the particle 
number densities: 

J^{t) = / dpn(p,i)= / dxn{x,t) . (5.45) 



The charge densities, the particle number densities as well as the total number of 
created particles are investigated in detail in Section 5.3. 



^^The momentum space particle number density n{p, t) is related to the one-particle distribution 
function J^{p, t) in the limit of a spatially homogeneous, time-dependent electric field E{t) according 
to Eq. (4.67): 

T{p,t) = '"l''(P>')-^vac(p)]+p|w(p,t)-Wvac(p)] ^ j;^ fn(p,t) . 
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5.2 Solution strategies 

The subsequent investigation of tlie Schwinger effect in space- and time-dependent 
electric fields E{x,t) is based on the equations of motion for the equal-time Wigner 
components Eq. (5.22) - (5.25), which are solved in order to calculate the charge 
densities, the particle number densities and the total number of created particles. 

As it is in general not possible to find an analytic solution for the equal-time 
Wigner components w{x,p,t), these calculations have to be based on numerical 
simulations. It has to be emphasized that the PDE system Eq. (5.22) - (5.25) is in 
fact linear and first order in the time derivative ^ as well as in the spatial derivative 

however, due to the appearance of the pseudo-differential operator: 

1 

A{x,p,t) = e r^diE[x + ii-l^,t)§^, (5.46) 

■'~2 

arbitrarily high momentum derivatives have to be taken into account. Accordingly, 
it is not possible to numerically solve this hyperbolic PDE system without further 
approximations or manipulations. 

Before actually discussing various solution strategies, it is convenient to trans- 
form the homogeneous PDE system Eq. (5.22) - (5.25) into an inhomogeneous one 
by defining modified equal-time Wigner components: 

m'"{x,p,t) = w{x,p,t) - Wvac(p) , (5.47) 

with: 

w"(x,p,tvac) =0 . (5.48) 

The equations of motion for the modified equal-time Wigner components w^{x,p,t) 
are accordingly given by: 

[g + A{x,p,t)]s- -2pp^ =A{x,p,t)^^, (5.49) 



[| + Zi(x,p,t)]vS =0 , (5.50) 

[|+Zi(x,p,t)]v^ +2mp^ = A{x,p,t)^ , (5.51) 

[§^+ A{x,p,t)]p'' +2ps" -2mv''=0 . (5.52) 



""^^It is assumed that the electric field E{x, t) can be expanded in a Taylor series with respect to 
the spatial coordinate. 
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5.2.1 Derivative expansion in p-space 



The expansion of the pseudo-differential operator Eq. (5.46) in a Taylor series with 
respect to the spatial coordinate reads: 



^-4"(2n+l)2n! 



E 



(2n). 



& 



2n+l 



dp 



2n+l 



(5.53) 



where E^^\x,t) denotes the k-tli derivative with respect to the spatial coordinate. 
Assuming that: 



E{x,t) 



dw''(x,p,t) E^^\x,t)d^w''{x,p,t) 



dp 



» 



24 



dp'^ 



(5.54) 



it is well justified to neglect the higher momentum derivatives so that A{x,p,t) 
might be approximated according to: 



d 

A{x,p,t)^eE{x,t)— . 



(5.55) 



Note that this expression becomes exact in the limit of a spatially homogeneous 
electric field. Additionally, due to the fact that the k-th derivative with respect to 
the spatial coordinate scales as: 



E^^\x,t) 



Ac ] ^ 

T 



(5.56) 



this approximation is assumed to be valid for weakly varying electric fields with the 
spatial variation scale being much larger than the Compton wavelength, A ^ Ac. 
On the other hand, it is clear that this approximation is expected to fail once 
spatial variation scales of the order of the Compton wavelength A ~ C(Ac) are 
under consideration. 

Accordingly, the equations of motion for the modified equal-time Wigner com- 
ponents Eq. (5.49) - (5.52) turn into a linear, first order hyperbolic PDE system in 
p-space within this approximation: 



[|+ei?(x,i)|]^s- 
[|+ei?(x,t)|]vS 
[|+ei?(x,f)|]v'' 



-eE{x,t)^ 




eE{x,t)l^]p^ 



gjwo +2mp^= eE{x,t)Jj^ 

+ 2ps'' -2771^"= 



(5.57) 
(5.58) 
(5.59) 
(5.60) 



with {x,p,t} S IR as well as w'"(x,p, tyac) = 0. 
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5.2.2 Full solution in y-space 

The range of validity of tlie derivative expansion is restricted as it is expected to 
fail for spatial variation scales of the order of the Compton wavelength A ~ 0{Xc)- 
Accordingly, a full solution of Eq. (5.49) - (5.52) is still needed in order to determine 
the actual range of validity of the derivative expansion as well as to check the quality 
of this approximate solution. Note that this full solution can surely be not obtained 
in p-space as arbitrarily high momentum derivatives occur when expanding the 
pseudo-differential operator A{x,p,t) for a general electric field E{x,t). 

Due to the fact that the momentum variable p appears linearly in the equations 
of motion for the modified equal-time Wigner components Eq. (5.49) - (5.52), it 
seems to be advantageous to transform these equations to conjugate y-space and 
solve them there. Along with the definition of the Fourier transformation: 



"^P e*PV(a;,p,t) , (5.61) 



(27r) 



momentum derivatives ^ transform into linear factors of y and vice versa: 



Note that: 



FT[pw''(x,p,i)] = -i^^-{x,y,t) , (5.62) 
FT[^w''(x,p,t)] = -iy^^x,y,t) . (5.63) 

w^{x,y,t) = [w''{x,-y,t)]* , (5.64) 



due to the fact that w'"{x,p,t) G R. Accordingly, one can restrict oneself to ?/ S 1R+. 
Additionally, the pseudo-differential operator Eq. (5.46) turns into a function of y 
upon Fourier transforming it:^'^ 



1 



FT [A{x,p, t)] =-iey f \ d^E{x + ^y,t) = -ieS{x, y, t) . (5.65) 

Finally, it is also necessary to Fourier transform the inhomogeneous part in order to 
formulate the PDE system Eq. (5.49) - (5.52) in y-space. This seems impossible on 
a first view as the Fourier transformation of Vvac(p) does not exist. Note, however, 
that the relevant quantity is in fact: 

FT [A{x,p, t)Wvac(p)] = ^^i^^ FT [|wvac(p)] , (5.66) 



^^The function £{x, y, t) has been introduced for later use. Note that the parameter integral over 
^ cannot be carried out analytically in general. 
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so that the remaining Fourier transformation can indeed be calculated in terms of 
modified Bessel functions of the second kind;^^ 

FT[|;:^] =-z^i^o(mM), (5.67) 

FT[|4y] = (5.68) 

Therefore, the equations of motion for the modified equal-time Wigner components 
Eq. (5.49) - (5.52) turn into a linear, first order hyperbolic PDE system in y-space: 

i'^£{x,y,t)Ko{my) , (5.69) 

, (5.70) 

^£{x,y,t)Ki{my) , (5.71) 

, (5.72) 

with y G 1R+ and {x,t} G IR as well as {x,y,tvax:) = 0. It has to be empha- 
sized that the solution of this PDE system is indeed a full one as compared to the 
approximate solution of the PDE system Eq. (5.57) - (5.60).^^ 

Note that the PDE systems which are obtained either from the Fourier transfor- 
mation to y-space or from a leading order derivative expansion in p-space are both 
linear and first order. Nonetheless, the numerical cost for actually solving them is 
quite different: As both of them are hyperbolic PDE systems they describe in fact 
propagation phenomena in {j;,p}-space and in {x, i/}-space. respectively. The speed 
of propagation can be deduced from the prefactors of the partial derivatives: 

y - space : ^ and 2^ , 

p- space: ^ and eE{x,t)-§^ . 

Considering the PDE system in y-space, there is in fact continued propagation in 
both x- and y-directions so that it is necessary to adopt an expanding grid in both 
directions for actually solving the PDE system numerically. For the PDE system 
in p-space, however, there is continued propagation only in x-direction whereas 
propagation in p-direction takes place only when the electric field E{x,t) is non- 
vanishing. Accordingly, an expanding grid is only needed in x-direction whereas one 
can stick to a fixed grid in p-direction. 

^*For details on Bessel functions see [113], chapter 9. 

^^An inverse Fourier transformation has to be performed in order to recover the modified equal- 
time Wigner components m" {x,p,t) in p-space. 



[|-ie^>,y,i)]i^S 

[| - ie£{x, y, t)] + f +2m = 

[| - ie£{x, y, t)] -2i-§- -2m = 
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5.2.3 Local density approximation 

The local density approximation is in fact no solution of the PDE system Eq. (5.49) 
- (5.52) but rather gives an approximation of the momentum space particle number 
density n[p,t) as well as the total number of created particles J^{t). This approxi- 
mation is based on the one-particle distribution function J-{p, t) which is obtained 
from solving the quantum Vlasov equation Eq. (4.36) - (4.38) in the presence of a 
spatially homogeneous, time-dependent electric field E{t). 

The idea is the following: Given that the spatial variation scale is much larger 
than the Compton wavelength, A ^ Ac, it is well justified to describe the Schwinger 
effect at any point Xi independently. Considering a separable space- and time- 
dependent electric field: 

E{x,t) = Eog{x)h{t) , (5.73) 

one can repeatedly solve the quantum Vlasov equation for a spatially homogeneous 
electric field E{t) with field strength E^gixi) at any point Xi. This yields the one- 
particle distribution functions J-{p,t;EQg{xi)), respectively, where the dependence 
on the field strength has been explicitly indicated. Correspondingly, one defines the 
particle number quasi-distribution in local density approximation: 

nioc(x,p, t) = 2J^{p, t; Eog{x)) . (5.74) 

Integrating over real space then gives the momentum space particle number density 
in local density approximation: 

nioc{p,t) = J j^nioc{x,p,t) , (5.75) 

whereas integration over the whole phase space yields the total number of created 
particles in local density approximation: 

Moc(i) = J dTni^,{x,p,t) . (5.76) 

It is clear that this approximation has a restricted range of validity as one has 
to assume that A ^ Xc so that the pair creation process at any point Xi can be 
considered as taking place in a spatially homogeneous, time-dependent electric field 
with field strength EQg{xi). Additionally, it has to be emphasized that no predictions 
on the charge densities can be made as there is local charge neutrality in the case 
of a spatially homogeneous, time-dependent electric field E{t). 
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5.3 Numerical results 

In this section I finally present the results of an ab initio simulation of the Schwinger 
effect in a simple space- and time-dependent electric field E{x, t) in QED^^^. Solving 
the equations of motion for the modified equal-time Wigner components allows 
to analyze the time evolution of the system and to investigate various observable 
quantities such as the charge densities, the particle number densities as well as the 
total number of created particles. This study should consequently shed some light 
on the effect of spatial and temporal variations on the pair creation process. 

As this is the first study of this kind, I restrict myself to a rather simple space- 
and time-dependent electric field: 

E{x,t) = Eog{x)hit) = Eoexp sech^ (^^^ , (5.77) 

where t determines the temporal extent of the pulse whereas A describes the charac- 
teristic length scale of the problem. The PDE system Eq. (5.57) - (5.60) originating 
from a leading order derivative expansion in p-space as well as the PDE system 
Eq. (5.69) — (5.72) derived from the Fourier transformation to y-space are then 
solved by means of a finite difference scheme with second order accuracy [125, 126]."'^® 
In order to actually solve the PDE system in t/-space later on, one additionally 
needs to calculate S{x,y,t) as defined in Eq. (5.65). The parameter integral of the 
spatial dependence over £^ is given by: 

'_.,..(fe±M).^.!lMMM), 



2 



SO that 



17 



£ix,y,t) = f-E,X [Erf (h^) + Erf (^)] sech^ (^) . (5.79) 

The parameters are chosen in such a way that the pulse length parameter takes the 
value T = ^ and the field strength parameter is given by e = 0.5. Accordingly, 
the Keldysh parameter 7 = 0.2 is still in the non-perturbative regime. The spatial 
variation scale A is considered for both A ^ Ac as well as for A ~ 0{Xc)- 



^^Details can be found in Appendix B. 

^^In order to solve the PDE system in j/-space efficiently, it is necessary to choose the spatial 
dependence in such a way that the parameter integral over ^ can be carried out analytically. 
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5.3.1 Time evolution 

This part is dedicated to the investigation of the time evolution of the marginal 
distributions for A = lOAc.-*^^ These observables are calculated by numerically 
solving the PDE system Eq. (5.69) - (5.72) in j/-space and subsequently performing 
an inverse Fourier transformation:"'^^ 



In Fig. 5.1, the time evolution of the momentum space marginal distributions is 
shown. These observables show several features which are partly known from the 
investigation of the Schwinger effect in spatially homogeneous, time-dependent elec- 
tric fields E{t): 

First of all, the pair creation process shows a non-trivial momentum dependence. 
Most notably, field excitations are created with momenta around p = and are then 
accelerated by the electric field to higher momenta. It has again to be emphasized 
that the interpretation of the particle number density n(p, t) as momentum distri- 
bution of real particles is only possible at asymptotic times t — > oo. At intermediate 
times, on the other hand, it can only be considered as mixture between real and 
virtual excitations. 

Moreover, the particle number density n(p, t) takes a fixed value once the electric 
field vanishes at asymptotic times i — > oo since there are no other forces acting than 
the electric.^'' Note that this behavior would change as soon as the self-induced 
electric field due to the pair creation process is taken into account. 

Additionally, one would naively expect that particles with charge e are accel- 
erated to positive momenta whereas particles with charge — e are accelerated to 
negative momenta. It is a peculiarity of the equal-time Wigner formalism, however, 
that both types of particles are accelerated to positive momenta. Consequently, a 
particle with charge — e and momentum p has to be interpreted as an antiparticle 
with physical momentum —p?^ This is further discussed when considering the time 
evolution of the real space marginal distributions in a moment. 

Finally, the charge density q{p,t) vanishes identically at all times. This can in 
fact be traced back to the choice of a symmetric spatial dependence g{x) = g{—x). 



It would also be possible to investigate the time evolution of the equal-time Wigner components. 

It is, however, much more enlightening to consider observable quantities instead of auxiliary ones. 
^®As the modified equal-time Wigner components are calculated on a finite grid in both x- and 

y-directions, one has to approximate the inverse Fourier transformation by a finite sum. 

^"Note that t = 3r is already very close to t ^ oo for the temporal dependence h{t) = sech^ {t)' 
^'^This peculiarity is well known from relativistic quantum mechanics, where a negative energy 

solution with momentum p is considered as an antiparticle with physical momentum —p as well. 
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Figure 5.1: Time evolution of the momentum space particle number density n{p,t) (blue) 
and charge density q{p,t) (purple). Note that another scale is used at t = 0. All other 
parameters are given in the text. 
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In Fig. 5.2, the time evolution of the real space marginal distributions is shown. 
Along with the discussion of the momentum space marginal distributions, a co- 
herent picture of the Schwinger effect in space- and time-dependent electric fields 
E{x,t) can be drawn: 

First of all, the investigation of the real space marginal distribution at early 
times t < shows that the pair creation process takes only place in the space region 
around x = where the electric field acts. It has to be emphasized that the spatial 
extent is actually determined by the spatial variation scale A.^^ 

Moreover, the particle number density n{x, t) is peaked around x = whereas the 
charge density q{x,t) vanishes there at early times t < 0. This can be interpreted 
as local charge neutrality in the center of the pulse whereas one observes charge 
separation at the edges of the pulse. The effect of charge separation is in fact driven 
by the electric field which accelerates excitations with charge e into the positive 
x-direction whereas excitations with charge — e are accelerated into the opposite 
direction. 

The effect of charge separation can be seen even better at later times t > 0. 
Owing to the acceleration in the electric field, one bunch of excitations with charge 
e propagates into the positive x-direction whereas another bunch of excitations 
with charge — e propagates into the opposite direction. These bunches can in fact 
be identified with particles and antiparticles, respectively, once the electric field 
vanishes either at asymptotic times t — > c» or beyond the region where the electric 
field acts. 

It has been mentioned when discussing the momentum space marginal distri- 
butions that both excitations with charge e and excitations with charge — e are 
accelerated to positive momenta against naive expectations. Nonetheless, they ac- 
tually move apart from each other in opposite directions even though both types 
of excitations do have positive momenta. This again motivates why a particle with 
charge — e and momentum p has rather to be interpreted as an antiparticle with 
physical momentum —p. 

Finally, the particle number density n(x,t) is symmetric with respect to reflec- 
tions at X = whereas the charge density q{x,t) is antisymmetric. This symmetry 
properties of the real space marginal distributions can again be traced back to the 
choice of a symmetric spatial dependence g{x) = g{—x). 



^^For the spatial dependence g{x) = exp( — ^^), the electric field has dropped to approximately 
1% of its maximum value at \x\ ~ 3A. 
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Figure 5.2: Time evolution of the real space particle number density n{x,t) (blue) and 
charge density q{x, t) (purple). Note that another scale is used at t = 0. All other parame- 
ters are given in the text. 
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5.3.2 Particle number density 

The time evolution of the marginal distributions has been considered in order to gain 
insight into the dynamics of the pair creation process in the presence of a space- and 
time-dependent electric field E{x,t). The investigation of the real space marginal 
distributions made in fact plain that bunches of particles and antiparticles move 
apart from each other in opposite directions. 

Note, however, that the dynamics of this motion is in fact governed by the 
momentum of the particles which is encoded in the momentum space marginal 
distributions. Due to the fact that the momentum space charge density q{p,t) 
vanishes identically at all times, it is indicated to focus on the momentum space 
particle number density n{p,t) in the following. 

5.3.2.1 Dependence on A 

This part is dedicated to the investigation of the dependence of the particle number 
density n{p, t) on the spatial variation scale A. To this end, the time evolution of 
n{p,t) for different values of A is investigated. Again, these results are obtained by 
numerically solving the PDE system Eq. (5.69) - (5.72) in y-space and performing 
an inverse Fourier transformation afterwards. For the sake of better comparability, 
it is in fact more convenient to investigate the reduced particle number density n{p, t) 
instead of n{p, t):'^^ 



In Fig. 5.3, the time evolution of fi{p, t) is displayed for different values of A, whereas 
the asymptotic value n{p, oo) is shown in more detail in Fig. 5.4. Disregarding the 
trivial scaling effect it turns out that n{p,t) shows various remarkable features as a 
function of A: 

First of all, the behavior of n{p,t) for different values of A is rather similar at 
early times t < 0, at least for A > 5Ac'; Field excitations are created with momenta 
around p = and are then accelerated by the electric field to higher momenta. 
Depending on the actual value of A, however, the shape of the asymptotic reduced 
particle number density n{p, oo) changes in various respects. 



Decreasing the spatial variation scale A results in a diminisliment of the total energy of the 
electromagnetic field £em{t) as well: 



Accordingly, there is less energy available and fewer particle-antiparticle pairs are created. This 
trivial effect of the decline of n{p, i) with decreasing A is in fact accounted for by considering n{p, i) 
instead of n(p, t). 




(5.81) 
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Figure 5.3: Time evolution of the reduced particle number density n{p, t) for various values 
of A: 100 Ac (blue), lOAc (purple), 5Ac (brown), Ac (yellow). Note that another scale is 
used at t = 0. All other parameters are given in the text. 
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Figure 5.4: Asymptotic reduced particle number density h{p, oo) for various values of A: 
lOOAc (blue), lOXc (purple), 8Xc (red), 6Ac (orange), 5Ac (brown), 4Ac (gray), 3Ac (dark 
green), 2Xc (light green), Ac (yellow). All other parameters are given in the text. 

Most notably, a decreasing value of A involves a shift of the peak momentum to 
smaller values: 



This behavior is understood in the following way: The field excitations are acceler- 
ated in such a way that the peak momentum takes a certain value pSo in the case 
of a spatially homogeneous, time-dependent electric field E{t). For a space- and 
time-dependent electric field E{x,t) = E{t)g{x), however, the value of the accel- 
eration by the electric field depends on the actual position x. Most notably, the 
field excitations are less accelerated compared to the spatially homogeneous case 
as \g{x)\ < 1. This finally results in a shift of the peak momentum jj^ to smaller 
values for decreasing values of A. 

Moreover, the shape of n{p, oo) becomes higher and narrower for decreasing 
values of A, at least for A > 4Ac- This is kind of a self- focussing effect which is 
caused by the spatial inhomogeneity. Excitations which are already created with 
high momenta are accelerated for a shorter period compared to excitations which 
are created with small momenta, which are then accelerated for a longer period. 
Consequently, the created particles are bunched into a smaller phase space volume. 

This behavior changes again for A < 4 Ac when n{p, oo) becomes in fact even 
narrower, however, the height of the peak momentum p^ decreases again. This is a 
direct consequence of the decreasing values of A as more and more excitations gain 
too less energy in order to finally turn into real particles. This behavior is most 
clearly seen for A = Xc when the energy content of the electric field is too small so 
that none of the field excitations turns into real particles at the end. 



Pai < Pxl Ai < A2 . 



(5.82) 
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5.3.2.2 Full solution vs. approximate solutions 

As it has become clear that n{p, t) shows a strong dependence on the spatial variation 
scale A especially for small values of A, it is about time to compare the full solution 
with the approximate ones. This comparison is actually done for two reasons: 

On the one hand, the comparison of n[y](p, t), which is obtained by solving the 
PDE system Eq. (5.69) - (5.72) in y-space, with n\oc{p,t) in local density approx- 
imation is suggested from a physical point of view: A deviation of ?T,jy](p,t) from 
flioc(p, t) indicates a change in the pair creation behavior in the sense that the pair 
creation process at any point Xi cannot be considered as taking place in a spatially 
homogeneous, time-dependent electric field E{t)g{xi) anymore. 

On the other hand, the comparison of n^y^{p^ t) with n[p](p, t), which is obtained 
by solving the PDE system Eq. (5.57) - (5.60) originating from a leading order 
derivative expansion, is rather important from a computational point of view as the 
numerical solution in y-space is much more expensive than the numerical solution 
in p-space. 

In Fig. 5.5, the full solution ^[y](p, oo) is compared with the leading order derivative 
expansion result n[p](p, oo) as well as with the local density approximation result 
fiioc(p, oo) for different values of A: 

Most notably, the difference between the various results is rather minor for 
A > lOOAc". Accordingly, the pair creation process at any point indeed 
be considered as taking place in a spatially homogeneous, time-dependent electric 
field E{t)g{xi) for large values of A. 

For decreasing values of A, however, the different results differ from each other: 
Most notably, the leading order derivative expansion result n[p] (p, oo) shows trun- 
cation artefacts for momenta around p = compared to the full solution n^y^[p, oo). 
Additionally, a shift of the peak momentum pf-^ < p^^^ occurs as discussed in detail 
previously. 

As expected, the leading order derivative expansion result fi^j{p, oo) becomes 
worse for decreasing values of A. However, it has not been anticipated that the 
leading order derivative expansion would fail for momenta around p = as a previous 
study had estimated that it would rather fail for large momenta [33].^'^ The leading 
order derivative expansion eventually breaks down once the truncation artefacts 
become even larger than the reduced particle number density itself. 

^''in fact, a different type of electric field had been considered in this investigation: 

E{x,t) = i;ocos(f )sech^(:^) . 
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Figure 5.5: Comparison of n[y] (p, cx)) (blue) with n[p](p, cx)) (purple) as well as nioc(p, oo) 
(red). Top: A = lOOAc- Middle: A = lOAc. Bottom: A = 5Ac. All other parameters are 
given in the text. 

5.3.2.3 Asymmetry of nyy^{p, oo) 

Even though ny^{p, oo) seems to be symmetric aromid the peak momentum p^*^, it 
turns out that there is in fact a small asymmetry which is displayed in Fig. 5.6 for 
two different values of A: 
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Figure 5.6: Asymmetry of the reduced particle number density around the peak momentum 
p^*^ with ^^[J,](^'^'' +p, oo) (blue) and nyy-^{p^^ —p,oo) (purple). Top: X = lOOXc- Bottom: 
X = lOAc. All other parameters are given in the text. 



First of all, this asymmetry seems to be a peculiarity of n^yj{p, oo) as it appears for 
all spatial variation scales A in such a way that: 

n[y] {'[f^ + oo) > n\y^ {p^^ - p, oo) with p > . (5.83) 

Additionally, one nicely observes again that the main peak becomes higher and nar- 
rower for decreasing values of A. This has already been discussed in detail previously. 

In order to better understand the origin of this asymmetry, it is convenient 
to consider the local density approximation for a moment. Due to the fact that 
fl[y](p, oo) and nioc(p, oo) nearly coincide for A > lOOAc, it should be possible to 
give an explanation already by considering the analytic expression for the asymptotic 
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Figure 5.7: Asymptotic one-particle distribution function !F{p, oo; EQg{xi)) for Xi = 
(blue), Xi = J (purple), Xi = ^ (red), Xi = ^ (orange) and Xi = A (brown). 



one-particle distribution function Eq. (4.145):^^ 

sinh [2ei?oT"+wo(p)-wi(p)] ) sinh [2e_BoT-aJo(p)+<.^l(p)] ) 



J='(p, oo; Eo) 



sinh \^7TTUJo{p)^ sinh {^7rTuJi{p)j 



with 



(5.84) 



^oip) = V "T'^ + b ~ 2eEQT]'^ and t^i(p) = \/ + . 



(5.85) 



Accordingly, the local density approximation result nioc(p, oo) might be approxi- 
mated by a Riemann sum: 



2 Ax 



N 



T{p, oo; -Bo) + 2 ^ .F(p, oo; Eog{iAx)) 



i=l 



(5.86) 



Note that each individual contribution J-{p, oo; E^gixi)) is in fact symmetric around 
its corresponding peak momentum: 



(5.87) 



Moreover, the peak momentum p'^{xi) of each individual contribution is shifted to 
a smaller value for increasing values of Xi as displayed in Fig. 5.7: 



for 



Xl > X2 ■ 



(5.88) 



^^Note that the asymptotic one-partice distribution function is now expressed in terms of the 
phase-space kinetic momentum p. Additionally, the dependence on Eo has been indicated exphcitly. 
Moreover, the factor 2 has been dropped since fermions are spinless particles in 1 + 1 dimensions. 
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"loc(P.O°) 




Performing the Riemann sum as displayed in Fig. 5.8, the peak momentum p^^^ is 
shifted to a smaller value compared to the peak momentum p^{0) of the one-particle 
distribution function J-{p,oo\Eq): 

pfoc<P?^'(0). (5.89) 

Accordingly, nioc(p, oo) becomes asymmetric around j^^^ as the individual contribu- 
tions J- {p, oo; EQg{xi)) are in fact asymmetric around p^^^^ as well. 

5.3.3 Total number of created particles 

I finally turn to the investigation of the number of created particles J^{t) . Actually, 
as the particle interpretation of the field excitations is only possible at asymptotic 
times t — > oo, I restrict myself to the asymptotic number of created particles A/"(oo). 

In Fig. 5.9, the full solution A/"[y](oo) is compared with the leading order deriva- 
tive expansion result A/'[p](oo) as well as with the local density approximation result 
A/ioc(oo) for different values of A. Note that it is again more convenient to consider 
the reduced number of created particles : 

moo) ^ ^ , (5.90) 



so that the trivial scaling effect with respect to A is disregarded. These results 
perfectly fit into the picture of the Schwinger effect in space- and time-dependent 
electric fields E{x,t) as concluded from the analysis of the marginal distributions: 
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Figure 5.9: Comparison of AT jy](oo) (blue) with A/'[p](cx)) (purple) and A/'io(.(cx)) (red) for 
different values of A in a lin-log plot. All other parameters are given in the text. 

First of all, A/'ioc(oo) shows an exact scaling behavior as function of A. This exact 
scaling exhibits itself as a straight line in the lin-log plot. The outcome of either 
method is then nearly the same for large values of A. This good agreement between 
the full solution with the approximate ones for large values of A has already been 
pointed out previously when discussing the results for the reduced particle number 
density n(p, oo) within the various schemes. 

For smaller values of A, however, the various solutions differ significantly. Most 
notably, both the full solution J7[y] (oo) and the leading order derivative expansion 
result A/'[p](oo) do not show the exact scaling behavior of A/'io(.(oo). Moreover, both 
of them predict a decrease in the total number of created particles compared to the 
local density approximation result. This behavior can be traced back to the fact 
that for a decreasing value of A more and more excitations gain too less energy in 
order to finally turn into real particles. 

Nonetheless, there is still a difference between the full solution Ar[y](oo) and the 
leading order derivative expansion result Ar[p](oo) for small values of A: Both of 
them show a significant decrease compared to ^^100(00), however, Ar[y](oo) actually 
drops to zero for A ~ Ac whereas Ar[p](cx)) does not. Actually, this is no surprise as 
it has already been shown that the leading order derivative expansion fails for small 
values of A. 

The sharp drop of Ar[y](oo) for small values of A is in fact in good agreement 
with previous studies of the Schwinger effect in space-dependent electric fields E{x) 
[22, 23, 109]. It has already been mentioned that the pair creation process terminates 
once the work done by the electric field over its spatial extent is smaller than twice 
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the electron mass: 



e J E{x)dx < 2m . (5.91) 
Specifically, for the space-dependent electric field: 

Eix) = Eoexp(^-^^ , (5.92) 

this condition reads: 

V2TTeEoX < 2m . (5.93) 

Accordingly, considering this space-dependent electric field E{x) for a field strength 
parameter e = 0.5, the pair creation process terminates for spatial variation scales: 

A < ij-Ac ~ 1.595Ac . (5.94) 
e V TT 

This estimate is in fact in good agreement with the nmnerical result for Ar[y](oo) in 
the presence of a space- and time-dependent electric field E{x, t) as well. 



Chapter 6 

Conclusions and outlook 



The main issue of this thesis was to investigate various aspects of the Schwinger 
effect in inhomogeneous electric fields. As a matter of fact, vacuum pair creation in 
the presence of simple field configurations such as static, sinusoidal or pulsed electric 
fields has been investigated since the early days of quantum mechanics, however, it 
was not until now that it becomes possible to study the Schwinger effect in realistic 
electric fields showing both temporal and spatial variations. This recent progress 
is due to a better theoretical understanding of the mechanism behind the pair cre- 
ation process as well as owing to the rapid development of computer technology 
which makes extensive numerical simulations feasible today. These advances in the 
theoretical description of the Schwinger effect are in fact urgently needed as a new 
generation of high-intensity laser systems such as the European XFEL or the Ex- 
treme Light Infrastructure are already in the starting blocks to shed some light on 
the strong-field regime of QED soon. 

The first part of this thesis was dedicated to the theoretical description of the 
Schwinger effect by means of quantum kinetic methods and briefly reviewed the 
Wigner formalism as well as the quantum Vlasov equation. In Section 4.2 it was 
eventually shown that the quantum Vlasov equation in its nowadays widely used 
form is in fact equivalent to the equal-time Wigner formalism in the presence of a 
spatially homogeneous, time-dependent electric field. This is quite remarkable as 
the quantum Vlasov equation is derived from canonical quantization whereas the 
equal-time Wigner formalism is formulated in phase space. As a consequence, it 
was even possible to calculate analytic expressions for the one-particle distribution 
function as well as for the equal-time Wigner components in the presence of both a 
static and a pulsed electric field in Section 4.3. 

It was a shortcoming of previous studies that the Schwinger effect has only 
been considered for simple time-dependent electric fields. This gap has been partly 
closed by the investigation of the pair creation process in the presence of a pulsed 
electric field with sub-cycle structure in Section 4.4. The corresponding results in- 
dicate that the momentum distribution of created particles is extremely sensitive to 
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the various laser parameters, resulting in a number of new observables such as the 
width of the momentum distribution or its oscillatory structure. As these character- 
istics might become crucial for the observation of the Schwinger effect in upcoming 
high-intensity laser experiments, future investigations should focus on the issue of 
pulse shaping in order to maximize the particle yield and provide clear signatures. 

A second drawback of previous investigations was that there have not been any rig- 
orous studies of the Schwinger effect in the presence of space- and time-dependent 
electric fields. In this respect. Section 5 provides substantial progress; Based on the 
equal-time Wigner formalism in QED]^_,_]^, an ab initio simulation of the Schwinger 
effect in such a type of electric field has been conducted for the first time, allowing 
for the calculation of the time evolution of various observable quantities such as 
the particle number density, the charge density or the number of created particles. 
These results provide in fact deeper insight into the pair creation process, however, 
they also raise further questions: 

The comparison between the full solution with a local density approximation 
showed that the momentum distribution strongly depends on the spatial variation 
scale A. Moreover, it turned out that the momentum distribution of particles is 
identical to that of antiparticles owing to the choice of a symmetric spatial variation 
g(x) = g{—x). In contrast to that, particles and antiparticles should behave differ- 
ently in the presence of an electric field which does not show this symmetry. It is, 
however, not clear a priori which consequences the choice of such an electric field 
would have. Accordingly, it would surely be worthwile to consider this question in 
future investigations as well. 

The comparison between the full solution with an approximate solution, which 
has been based upon a leading-order derivative expansion of the occurring pseudo- 
differential operator, was another issue of investigation. This study in fact showed 
that the approximate solution deviates significantly from the full solution already 
at comparatively large spatial variation scales A. Quite surprisingly, this deviation 
occurs at small kinetic momentum values and leads to a break down of the leading- 
order derivative expansion for spatial variation scales A which are still somewhat 
beyond the Compton wavelength A^. 

Even though this first investigation of the Schwinger effect in a space- and time- 
dependent electric field in QED]^_|_j^ represents substantial progress, one has to admit 
that this study has to be viewed only as the first step towards a full description of 
the pair creation process in upcoming high-intensity laser experiments. As a matter 
of fact, there are still various issues which should be tackled in future investigations: 
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First of all, the space- and time-dependent electric field which has been under 
consideration cannot be considered as realistic in the sense of representing a field 
configuration to be realized at high-intensity laser facilities. As a next step one could 
in fact try to consider a more realistic time dependence such as a pulsed electric 
field with sub-cycle structure again. 

Secondly, the numerical simulation which has been performed was based on the 
equal-time Wigner formalism in QED]^_,_]^. Accordingly, the notion of a magnetic field 
did not even exist and Dirac fermions had to be considered as spinless particles. 
In order to avoid these shortcomings one should rather consider the equal-time 
Wigner formalism in QED3^]^ in future investigations. However, as it does not seem 
feasible to solve the equations of motion for the equal-time Wigner components for 
an arbitrary field configuration for computational reasons in the near future, one 
should restrict oneself to highly symmetric configurations for now. 

Finally, as soon as it becomes doable to perform numerical simulations of the 
Schwinger effect in the presence of space- and time-dependent electromagnetic fields 
in 3+1 dimensions, one should also reconsider the backreaction issue. As a start one 
could in fact be satisfied with a mean field description, however, the ultimate goal 
should be to consistently describe the pair creation process beyond the mean field 
level by taking into account photon corrections to the background electromagnetic 
field. 



Appendix A 

Specific calculations 



This appendix is devoted to giving the details of some calculations. 

A.l Covariant Wigner formalism: Equation of motion 

The starting point for the derivation of the equation of motion for the covariant 
Wigner operator W(x, p) is the underlying gauge-invariant density matrix C(A;x,y): 

C(A; X, y) = U{k; x, y) [^{x^ = x - *(x2 = x + |)] . (A.l) 

Note that the derivatives with respect to xi and X2, respectively, can be expressed 
as: 

d;' = \dl-dl and (f^ = lcf^ + dl. (A.2) 

Taking the derivative of C(A; x, y) with respect to X2 will eventually result in the equa- 
tion of motion for the Wigner operator Eq. (3.10). Note that taking the derivative 
with respect to Xi instead would result in the adjoint equation of motion Eq. (3.11): 

[I] m 

7'^a;i^C(A;x,y) =d^lJ{A^^t' ['l'(xi), ^(xa)] +iY(A;x,y)7^ [^{x^),d;^^{x2)] . 

(A.3) 

In order to calculate part [I], one uses Eq. (A.2) and performs the derivatives with 
respect to x and y, respectively: 

1 

d^^l({A;x,y)=ie f {{^ + ^)y-d^A,{x + ^y) + A^{x + ^y)}l({A;x,y) . (A.4) 

Using the definition of the field strength tensor F^i^{x) = d^Ai^{x) — d,^Af^{x) along 
with the identity g|^^(x -|- ^y) = y'^d^A^{x + ^y), the integral can be rewritten as: 

1 

/'ide{(i+6^^,.(x + ey)2/" + ^ [(i+C)A^(x + ey)]}. (A.5) 
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Accordingly, part [I] is given by: 



U{A;x,y). (A.6) 



Part [II], on the other hand, can be calculated immediately by taking into account 
the Dirac equation Eq. (3.2): 

t'' [f'(xi),9;:,^*(x2)] = -im [*(xi),^'(x2)] -iejf'A^ix2) [f'(xi), ^(xa)] . (A.7) 

Adding up the both parts [ I ] and [ II ] , one obtains the equation of motion for the 
gauge- invariant density matrix:^ 



i^,,(x, y) + n^(x, y) 7^C(A; x, y) = -imC{A: x, y) 



(A., 



with the operators D^{x,y) and n(x, y) being given by: 

1 

D^(x,y) = d^^-ie [\d^F^,{x + ^y)y\ (A.9) 

(A.IO) 



n^(x, y) = dl-ie I daF,,,(x + iy)y- . 

•'~2 

In order to switch to momentum space, one performs a Fourier transformation with 
respect to y. Accordingly, derivatives with respect to y become linear factors in p 
and vice versa: 

dl -iPii. and -id^ . (A. 11) 



This immediately results in the equation of motion for the covariant Wigner operator 
Eq. (3.10):^ 

[\D^{x, p) - in^(x, p)] 7^W(4)(x, p) = -imW^^\x, p) . (A.12) 



^Taking the derivative of C(x,y) with respect to xi instead of X2 yields; 

[ir»^(x, y) - fl^(x, y)] C(A; x, 7)7" = imC(A; x, y) . 

^Taking the Fourier transformation of the adjoint equation yields the equation of motion 
Eq. (3.11): 

[ii?,,(x,p) + in^(x,p)] W(^)(x,p)7'' = imW('''(x,p) . 
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A. 2 Equal-time Wigner formalism: Equation of motion 

An energy average is performed in order to switch from the covariant Wigner for- 
malism to the equal-time Wigner formalism. Considering the pseudo-differential 
operators Eq. (3.12) - (3.13):^ 

1 

2 



Do{x, p) =d^, + e I deE(x - i^dp) ■ Vp , (A. 13) 

1 

D(x, p) = Vx + e dC [E(x - i^dp)dp, + B(x - i^dp) x Vp] , (A.14) 

1 

no(x, p) = po + ie l\ deeE(x - i^dp) ■ Vp , (A. 15) 

1 

n(x, p) = p-ie f\ d^^ [E(x - i^dp)dp, + B(x - i^dp) x Vp] , (A.16) 

one assumes that the electromagnetic field can be expanded in a Taylor series with 
respect to the temporal coordinate: 

oo 

F^,(x + i^Vp, xo - i^dp,) = h.^l^i^ + ^?Vp, XQ){-iidp,r ■ (A.17) 

n=0 

(n) 

Here, F^i,^ denotes the n-th derivative with respect to the temporal coordinate. 
Upon taking the energy average of any of these operators acting on the covariant 
Wigner components W(x, p), one finds: 

/ Do{x, p)W(x, p) = A(x, p, i)w(x, p, t) , (A.18) 



J ^ D(x, p)W(x, p) = D(x, p, i)w(x, p, t) 



(A.19) 



no(x, p)W(x, p) = nt(x, p, i)w(x, p, t) + wW (x, p, t) , (A.20) 

/-^n(x,p)W(x,p) = n(x,p,t)w(x,p,i) , (A.21) 
J (27r) 

given that the covariant Wigner components and all its derivatives with respect to po 
vanish for \pq\ — > oo. The transport equations Eq. (3.36) - (3.43) and the constraint 
equations Eq. (3.48) - (3.55) are accordingly found by taking the energy average of 
the equations of motion for the covariant Wigner components Eq. (3.22) - (3.31). 



^Note that the operator D(x, p) corresponds to the covariant components _Di(x, p)ei whereas 
n(x, p) is related to the contravariant components ir(x, p)ej. 
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A. 3 Hierarchy truncation 



In order to derive the equations of motion for the first energy moments w^^^ (x, p, t), 
one first multiplies the equations of motion for the covariant Wigner components 
Eq. (3.22) - (3.31) by a factor pQ and subsequently takes the energy average; 



2U ■ t^^i -2nW ■ iti 



2n t 



[1] 



2U 



[1] 



t2 



Dt ' +d\^^ vo +D ■ vt^l +D[i1 ■ V 

Dt a|)^' +d\^^ ao +D ■ a^^l +d[^1 ■ a 

A vl^l +^1^^ ^ ^[1] _^2n X al^l +2n[^I X a = 

Dt al^l +1?^^ a +D aj,^' +D[i1 ao +2n x v^^] +2n[il x v = 

Dt t'^' ^ ti +D X t^^' +d[^i X t2 +2U ffii^i +2n[^i s - 

Dt t^2^ +d\^^ t2 -D X t^^' -DI^I X iti -2npW -2n[^l p = 





-2m a, 



[1] 





2m 

-2m t^^^ 


2m 





(A.22) 
(A.23) 
(A.24) 
(A.25) 
(A.26) 
(A.27) 
(A.28) 
(A.29) 



Here, the new operators D\^\x,p,t), Y)^^\x,p,t) and n[^](x,p,t) are given by: 
1 

DW(x,p,t) =ie daEW(x + ieVp,t)-Vp, (A.30) 

1 

DW(x,p,i) =ie de[CB«(x + ieVp,t) X Vp+ieE(x + ieVp,t)] , (A.31) 

1 

nW(x,p,i) = e d^[e'BW(x + ieVp,t) X Vp + ieE(x + i^Vp,t)] . (A.32) 



One can then express w[^l(x, p, t) in terms of w(x, p,t) by using the constraint 
equations Eq. (3.48) - (3.55). Various commutators are needed in order to simplify 
the resulting expressions when pulling through the various operators. Additionally, 
the homogeneous Maxwell's equations: 



Vx-B(x,i)=0 and b(^)(x, i) + Vx x E(x, t) = 



(A.33) 



as well as the representation of the electromagnetic field at shifted arguments are 
used:^ 



Pt'^U + i^Vp, t) = e'^^F^"'{x, t) = Y , F'^'i^, t 



(A.34) 



n=0 



The triangle operator is a shorthand notation for A = Vx • Vp 
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The relevant commutators including Dt{x,p,t) are given by: 

[A,nt] =i?W , (A.35) 
1 

[Dt , D,] = - e /"'^ A Ei(x + i^Vp.t) , (A.36) 

•'~2 

[A,n,] = Vx,int-Dfi . (A.37) 

The commutators including Ilt{x,p,t) read: 

[Ut,-Di] = - Vx,.nt , (A.38) 

1 

[n* ,n,] = - nfi + e f\ d([2i^ - e^W^ + ^^Vp,*) . (a.39) 

Finally, the commutators between D(x, p, t) and r[(x, p, t) are given by:^ 

1 

[D,,D,-] = e /"'^dCe.jfcABfc(x + ieVp,0 , (A.40) 

•'~2 
1 

[n,: , D,-] =-e f\ d^e^jk [1 + i?A] B;t(x + i^\/p,t) , (A.41) 

1 

[Ui,Uj] = e d^e,jk [2i^ - fA] Bk{^ + ^^\/p,t) . (A.42) 



As a matter of fact, one finds: 



1 



e dC[2ie-e2A + iA]F^'^(x + ieVp,t) = 0, (A.43) 

which can be easily checked by plugging in the Taylor expansion Eq. (A. 34) and 
calculating order by order the coefficients of A". Accordingly, the following identities 
hold: 

nf 1 + [Ut , n,:] - i [A , D,] = , (A.44) 
[ni,n,] +i[D,,D,-] = . (A.45) 



^In fact, the actual value of the commutator [11;, Dj] is only needed insofar as one needs to 
know: 

[n,,D,]=0 and [n,, D^] = [D,, n^] . 
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In order to simplify notation afterwards, I recall the equations of motion for the 
equal-time Wigner components Eq. (3.36) - (3.43) at this point:^ 



As 




- 2n ■ ti 




, 


(A.46) 


Ap 




+ 2n ■ t2 


+ 2mao = 


, 


(A.47) 


A vo 


+ D ■ V 






, 


(A.48) 


Dt ao 


+ D ■ a 




— 2mp = 


, 


(A.49) 


Av 


+ Dvo 


+ 2n X a 


+ 2mti = 


, 


(A.50) 


Dt a 


+ D ao 


+ 2n X V 




, 


(A.51) 


Ati 


+ D X t2 


+ 2nfB 


— 2mv = 


, 


(A.52) 


A t2 


- D X ti 


- 2np 




. 


(A.53) 



Taking into account all the identities Eq. (A. 35) - (A. 45) as well as the equations of 
motions for the equal-time Wigner components Eq. (A.46) - (A.53), the equations 
of motion for the first energy moments can be written as: 





Ihs 


Eq. 


{AA6)] 


+ iD ■ [Ihs Eq. (A.52)] 


= 0, 


(A 


54) 


-Ut 


Ihs 


Eq. 


(A.47)] 


- iD ■ [Ihs Eq. (A.53)] 


= 0, 


(A 


55) 


-lit 


Ihs 


Eq. 


(A.48)] 


+ n- [Ihs Eq. (A.50)] 


= 0, 


(A 


56) 


-Hi 


Ihs 


Eq. 


(A.49)] 


+ n- [Ihs Eq. (A.51)] 


= 0, 


(A 


57) 


-Hi 


Ihs 


Eq. 


(A.50)] 


- iD X [Ihs Eq. (A.51)] + n [Ihs Eq. (A.48)] 


= 0, 


(A 


58) 


-Hi 


Ihs 


Eq. 


(A.51)] 


- iD X [Ihs Eq. (A.50)] + U [Ihs Eq. (A.49)] 


= 0, 


(A 


59) 


-lit 


Ihs 


Eq. 


(A.52)] 


- iD [Ihs Eq. (A.46)] + n X [Ihs Eq. (A.53)] 


= 0, 


(A 


60) 


-Ut 


Ihs 


Eq. 


(A.53)] 


+ iD [Ihs Eq. (A.47)] - n X [Ihs Eq. (A.52)] 


= 0. 


(A 


61) 



Accordingly, it turns out that all these equations are trivially fulfilled once the 
equal-time Wigner components obey the equations of motion Eq. (A.46) - (A.53). 
Consequently, one does not have to solve the infinite hierarchy but it suffices to use 
the constraint equations in order to calculate all higher energy moments.''' 

®Note that the mass terms have been brought to the left hand side (denoted as 'Ihs' in the 
foUowing) compared to Eq. (3.36) - (3.43). 

'This calculation for first energy moments w'^'(x, p,f) served only as demonstration. As a 
matter of fact, it has been proven generally that the equations of motion for all the higher energy 
moments w'"'(x, p, t) are trivially fulfilled as well [93]. 
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A. 4 Vacuum value of the equal-time Wigner function 

In order to calculate the vacuum value of the equal-time Wigner function one canon- 
ically quantizes the free Dirac field as discussed in Section 4.1. Accordingly, one has; 

*(x,t) = J ^e^i-J] [uM,t)ar{q)+Vr{-q,t)bl{-q)] , (A.62) 

with:^ 

Wr(q,t) =[-/''io{q)-'yq + m]gitKfl,t)Rr, (A.63) 

Vri-q,t) = [-70cj(q) - 7-q + m]<7y(q,i)i?r • (A.64) 

The properly normalized positive and negative energy plane wave solutions gi^} (q, t) 
are given by: 

aiti (q, t) = ^ • (A.65) 

The spinors ^^(q, i) and ^^(q, i) obey the following completeness relations:^ 

Vwr(q,t)ur(q,t) = — — -[7°w(q) -7-q + m] , (A.66) 

Vz;r(q,i)t;r(q,t) = — — [7°a;(q) - 7 ■ q - m] . (A. 67) 

^ 2w(q) 

The equal-time Wigner function for the free Dirac field is defined as: 

Ma^(x,P,t) = IJ d'ye-^^-y{n\ [t(x- f ,t),*(x+ f |Q) . (A.68) 

Plugging in the explicit expressions Eq. (A.62) as well as using the non-vanishing 
anticommutators : 

{a,(q),4(qO} = {6,(q) , 6t(q0} = (27r)35,,5(q - q') , (A.69) 



*For simplicity, the same spinor representation as the one chosen in Section 4.1 is used. Note, 
however, that any other complete set of spinors would do the job equally well. 

^Note that the common factor 2o}(q) non-standard but a direct consequence of the chosen 
spinor normalization: 

uj.(q, t)us(q, t) = (q, t)vs{q, t) = 6ra ■ 
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one finds: 



I 2 

Wiac(x,P,i) = n XI K(-P,i)i^r-(-P,i) - Ur{p,t)Urip,t)] . 



(A.70) 



r=l 



Taking advantage of the completeness relations Eq. (A. 66) - (A. 67), this finally 
yields: 



m . 



(A.71) 



Accordingly, the vacuum value of the equal-time Wigner function is spatially homo- 
geneous and time-independent as it should be. Moreover, it turns out that sivac(p) 
and Vvac(p) are the only non- vanishing equal-time Wigner components: 



^vac (p) 



2m 
u;(p) 



and Vvac(p) 



2p 

uj{p) 



(A.72) 



A. 5 Basis expansion 



In order to expand ^(q, i) in the basis {ei,...,eio} as indicated in Eq. (4.57), it is 
convenient to choose the first basis vector in such a way that £^i(q, tyac) = 1 is the 
only non-vanishing coefficient in the vacuum: 



(Bi(q, t) 



a;(q,i) 



m 



7v{q,t) 


V / 



(A.73) 



Acting with M(q, t) as well as with the total time derivative on this first basis vector 
(Bi(q, t), one finds: 



(q,t)ei(q,i) = and <Bi(q,t) = -Q(q, i)<B2(q, i) , 



with the second basis vector e2(q, t) being defined as: 



(A.74) 



(E2(q,i) 



1 



rmrs{q3,t) 



7r(q,i)7r3(g3,i) - uj^{q,t)e3 





w(q,i)e_L 

Repeating this procedure with the second basis vector (B2{ci,t), one finds: 



(A.75) 



M(q,t)e2(q,i) = 2w(q,t)e3(q,t) and ii2iq,t) = Q{q,t)<Bi{q,t) , (A.76) 
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with the third basis vector (B3(q, being given by: 



<E3(q,i) 








Tv{q,t) X eg 
\ -meg y 



(A.77) 



Continuing with the third basis vector (B3(q, t), it turns out that the system closes 



M{q,t)<B3{q,t) = -2uj{q,t)<B2{q,t) and (E2(q,i)=0. 



(A.78) 



As £^i(q, tvac) = 1 is the only non- vanishing coefficient in the vacuum, one does in 
fact not need to calculate the remaining basis vectors. Accordingly, the system is 
fully described by means of the first three expansion coefficients Si{q,t), £2{^.,t) 
and £:3(q, t). 

A. 6 Series reversion 



One needs to invert the variable transformation in order to calculate the WKB 
instanton action 2<Sinst.(Q)' 



2= '^■''^(vf^^Vl' 



c.c. . 



(A.79) 



Due to the fact that the series expansion of the error function Erf(x + iy) is given 
by [113]: 



Erf (x + iy) 



2 

77^ ^ 



n!(2n + 1) 



iy) 



2n+l 



n=0 
^ oo 2n+l 

Hi ^ ^ n 

n=0 m=0 



(-1)" (2n + 1 
!(2n + l)V ra 



Jln+\-rai -xm 



iw)'' 



(A.80) 



the series expansion of T can be written as: 



e 2 



EE 

n=0 m=0 

1 + ^2 



-1 



2"n!(2n + l)T 



2n— 2m 



2n + 1 
2m 



2n+l-2m 



6t^ 



3 + + a' 



120r4 



t^+0(t^) . 



(A.81) 



As a matter of fact, the basis vectors (Ei(q, t) with i = {1, 2, 3} form an orthonormahzed set. 
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Given the series expansion of any function y{x) with y{0) = as well as y'{0) 7^ 0: 

00 

y{x) = Cnx"^ with Ci ^ , (A.82) 

n=l 

the inverse series expansion can be calculated term by term: 

oo 

x(2/) = 5]d„y" , (A.83) 

n=l 

with the first few coefficients being determined via [113]: 





= 0, 


cfd2 


= -C2 , 


4d3 


= - C1C3 , 




= 5C1C2C3 — C1C4 — , 


Cid5 


= 6C1C2C4 + 3CxC3 + 14C2 — Cj^Cs — 2IC1C2C3 



Accordingly, the inverse series t{T) can be calculated term by term by taking into 
account the series expansion Eq. (A. 81), 

m = T + 1±^!t3 + ^i^^^^T^ + OiT^) . (A.84) 



Appendix B 

Modified Lax-Wendroff scheme 



In order to solve the PDE systems Eq. (5.57) - (5.60) and Eq. (5.69) - (5.72), 
respectively, a modified Lex-Wendroff scheme for the solution of hyperbolic PDE 
systems with source terms has been used [125, 126]. A linear, first order PDE 
system in two dimensions reads in general:^ 

d d d 

— w(a, /?, t) = ^w(a, /3, t) + ^w(a, /?, t) + s(a, /3, i; w) , (B.l) 



with w(a, /3,i) being a n-component vector of field variables and s(a, /3,i;w) de- 
noting the n-component source vector: 



s(a, /?, t; w) = Cw(a, /3, t) + s(a, /3, t) . 



(B.2) 



The PDE system Eq. (B.l) is in fact hyperbolic given that all linear combinations 
of the n X n matrices A^^j and Bj^j are diagonalizable with real eigenvalues.^ The 
domain is then divided into rectangular cells in the a-/3 plane in order to actually 
solve this PDE system by means of a finite difference scheme. Moreover, it is 
convenient to introduce the following abbreviation for the vector of field variables 
at a grid point (a^, (3j) at a given time t„: 



(B.3) 



^The a- and /3-direction correspond to the x- and p-direction for Eq. (5.57) - (5.60) as well as 
to the X- and the y-direction for Eq. (5.69) - (5.72), respectively. 

^It can be checked easily that the PDE system Eq. (5.57) - (5.60) with: 



A. 



( 



V 



\ 

-1 



/ 



and B 



b] 



/ -eE{x,t) 



V 





-eE{x,t) 








-eE{x,t) 




as well as the PDE system Eq. (5.69) - (5.72) with: 



A, 



( 




V 



\ 

-1 



/ 



and -B[j,] = 



/ 





V 2i 



-2^ \ 



I 



\ 



-eE{x,t) J 



fulfill this hyperbolicity condition. 



110 



The numerical solution at a subsequent time tn+i = tn + At can for instance be 
found by adopting an explicit two-step scheme for discretizing the PDE system.^ 
Introducing the abbreviations: 

= ^■,;[.]wi:, , (B.4) 

s^, = a;:,wi;,. + s^, , (B.6) 
the first step utilizes the data w"" in order to calculate: 
w 

"'i+l/2,j + l/2 



(B.7) 





1 

4 






+ <+lj 




+ 


4 






+ ^r+i,j 


— a" 


+ 


4 


[ t>r+i j+i + 




— h" 


- b" 


+ 


At 
8 








+ s"- 



with: 



A[„i = ^ and Xtm = . (B. 



Aa ^nPj ~ A/3 • 

The second step, on the other hand, utilizes the intermediate data w"+^/^ in order 
to calculate:'* 

^n+l _ n , \c^i n+l/2 _ n+1/2 , a"+l/2 _ n+1/2 i 

ij 2 L'*i+l/2J + l/2 '*i-l/2j + l/2 "'"'*i+l/2j-l/2 '^i-l/2j-l/2 J 

A[^r, n+1/2 ,u"+l/2 _H"+l/2 _u"+l/2 1 

2 L i+l/2,j + l/2 ^ i-l/2j + l/2 i+l/2,j-l/2 i-l/2j-l/2j 



+ l/2,j + l/2 ' i-l/2,j + l/2 ' i+l/2,j-l/2 ' i-l/2,j-l/2 J ' 

(B.9) 



This two-step scheme is in fact consistent and of second-order accuracy. The stability 
condition reads: 

cA < 1 with A = max (A[q], A[^]) , (B.IO) 

with c denoting the absolute maximum of the propagation speeds which in turn is 
calculated by determining the absolute maximum of the eigenvalues of: 

A[a] cos((/<) + -B[^] sin((/)) with l^l < vr . (B.H) 



^The initial value problem is set up by choosing appropriate initial conditions at a certain initial 
time to- 

^Note that it is necessary to choose appropriate numerical boundary conditions at the edge of 
the domain in order to solve the initial value problem on a finite grid (ai,/3j). 
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